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MASSACHUSETTS 


EXECUTIVE  SUMMARY 


Introduction 

The  application  of  adaptive  optics  technology  to  laser-beam  control  has  been  under  investigation 
within  the  military  community  since  the  early  1970s.  This  work  has  encompassed  a  wide  range  of 
theoretical  investigations  as  well  as  an  impressive  suite  of  laboratory  and  field  experiments.  During  this 
period  a  number  of  real-time  systems  incorporating  several  hundred  degrees  of  freedom  have  been  built 
and  tested,  and  the  technology  to  extend  this  number  by  at  least  an  order  of  magnitude  is  now  under 
development. 

To  exploit  adaptive-optics  concepts  in  the  construction  of  astronomical  telescopes,  the  budgetary 
restrictions  under  which  this  field  of  research  operates  must  be  considered  and  improvements  in  system 
performance  must  be  carefully  balanced  against  equipment  complexity,  reliability,  and  cost.  To  accom¬ 
plish  this  design  trade-off,  a  set  of  concise  analytic  models  describing  the  essential  properties  of  a  laser 
guide-star  phase-conjugation  system  has  been  assembled.  With  the  aid  of  these  models,  an  optimization 
strategy  for  mating  adaptive  optics  to  a  4-m-class  optical  telescope  has  been  developed;  research  shows 
that  such  a  system  might  be  expected  to  improve  the  effective  atmospheric  seeing  conditions  by  nearly 
a  factor  of  10  within  the  isoplanatic  patch  of  the  turbulence  probe. 

The  Laser  Guide-Star  Concept 

Few  celestial  objects  are  bright  enough  to  permit  accurate  wavefront  measurements  by  a  phase 
sensor.  A  solution  to  this  problem  was  first  discovered  in  the  early  1980s  by  investigators  in  the  high- 
energy  laser-beam-control  community.  The  technique  involves  the  use  of  a  pulsed  laser  beam  that  is 
projected  along  the  path  to  the  target  object;  backscatter  generated  by  molecules  and  particulates  in  the 
atmosphere  provides  a  synthetic  probe,  which  can  be  used  to  obtain  reliable  turbulence  measurements. 
The  quality  of  the  data  provided  by  this  type  of  source  is  strongly  affected  by  the  altitude  of  the  beacon 
and  the  size  of  the  collection  aperture.  Residual  phase  errors  due  to  focal  anisoplanatism  (which  is  related 
to  the  difference  in  range  between  the  synthetic  beacon  and  the  target  object)  and  offset  anisoplanatism 
(which  arises  from  the  need  to  obtain  tracking  information  from  a  bright  neighboring  star)  represent 
driving  factors  in  the  overall  performance  of  the  adaptive-optics  system. 

An  Analytical  Treatment  of  Atmospheric  Turbulence 

To  a  large  extent,  the  success  of  this  investigation  can  be  attributed  to  the  extensive  use  of  a  new 
analytic  model  that  accurately  describes  the  principal  anisoplanatic  effects.  This  formulation  incorporates 
the  Rytov  approximation  (appropriate  for  slowly  varying  turbulence  and  low  scintillation)  in  the  devel¬ 
opment  of  integral  expressions  for  the  aperture-average  phase  error.  In  the  most  general  case,  the  solution 
requires  a  six-fold  integration,  but  four  of  these  integrals  are  easily  performed  using  standard  mathemati¬ 
cal  techniques.  Of  the  two  remaining  integrals,  which  involve  the  transverse  spatial  spectrum  and  the 
propagation  path,  the  first  is  performed  using  an  innovative  approach  that  applies  Mellin-transform 
theory.  For  most  problems  of  interest,  the  final  solution  can  be  expressed  as  a  rapidly  converging  series 
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expansion  that  incorporates  the  turbulence  moments.  A  brief  description  of  the  approach  and  derivations 
of  the  most  important  phase-error  formulas  can  be  found  in  Appendix  A. 

An  Optimized  Design  for  Guide-Star  Systems 

For  an  adaptive-optics  system  to  be  useful  in  astronomical  investigations,  a  significant  improve¬ 
ment  in  resolution  (in  excess  of  a  factor  of  2)  must  be  achieved  over  a  large  fraction  of  the  nighttime 
sky.  A  series  of  computer  simulations  in  which  the  quality  of  the  phase  correction  was  varied  has 
determined  that  the  high-spatial-frequency  (figure)  component  of  the  phase  error  can  be  no  more  than  3 
rad2  if  this  goal  is  to  be  achieved.  If  this  criterion  is  satisfied,  the  effective  resolution  of  a  long-exposure 
image  will  be  determined  by  tracking  jitter,  which  in  most  cases  is  dominated  by  anisoplanatic  errors 
related  to  the  angular  offset  between  the  target  object  and  the  tracking  star.  Therefore,  resolution  and  sky 
coverage  are  intimately  related;  an  improvement  in  one  of  these  factors  can  only  be  achieved  to  the 
detriment  of  the  other. 

Figure  ES-1  illustrates  the  relationship  between  achievable  resolution  and  fractional  sky  coverage. 
Although  it  is  theoretically  possible  to  build  an  adaptive-optics  system  that  provides  a  tenfold  improve¬ 
ment  in  the  resolution  of  a  4-m  telescope  operating  in  the  visible,  the  probability  of  finding  a  suitable 
tracking  star  to  maintain  the  requisite  tilt  jitter  is  predicted  to  be  only  1  percent.  On  the  other  hand,  a 
factor-of-4  improvement  can  be  attained  with  near-unity  probability.  As  discussed  in  Section  4,  such  a 
system  would  require  approximately  300  actuators  operating  at  a  correction  bandwidth  of  20  Hz  and  a 
beacon  laser  capable  of  maintaining  an  average  power  level  of  50  W. 


SYSTEM  RESOLUTION  (arc  sec) 


Figure  ES-I.  The  fraction  of  the  sky  over  which  a  given  resolution  can  be  attained  for  a  laser  guide-star  compen¬ 
sation  system.  Effects  of  residual  image  motion  due  to  finite  tracker  bandwidth,  limited  signal,  and  offset  anisoplanatism 
have  been  included.  The  uncorrected  seeing  at  0.55  pm  is  0.85  arc  sec  for  the  atmospheric  conditions  used  in  this 
study. 
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1.  INTRODUCTION 


The  application  of  adaptive  optics  to  laser-beam  control  problems  has  been  under  investigation 
within  the  military  community  since  the  early  1970s.  This  work  has  encompassed  a  wide  range  of 
theoretical  investigations  as  well  as  an  impressive  suite  of  laboratory  and  field  experiments.  During  this 
period  a  number  of  real-time  systems  incorporating  several  hundred  degrees  of  freedom  have  been  built 
and  tested,  and  the  technology  to  extend  this  number  by  at  least  an  order  of  magnitude  is  now  under 
development. 

The  designers  of  laser-propagation  systems  are  frequently  concerned  with  uncooperative-target 
engagements  in  which  the  target  return  is  either  too  weak  or  improperly  positioned  to  be  used  effectively 
as  a  turbulence  probe.  A  potential  solution,  which  involves  the  use  of  backscatter  from  laser  beams 
focused  in  the  atmosphere,  was  first  proposed  to  the  defense  community  by  Feinleib  in  1982  [1]  and  later 
independently  published  by  Foy  and  Labeyrie  in  1985  [2],  Since  that  time,  a  number  of  alternative 
schemes  have  been  considered  to  achieve  near-difffaction-limited  beam  profiles;  these  schemes  include 
arrays  of  low-altitude  sources  and  high-altitude  beacons  formed  within  the  earth’s  sodium  layer. 

This  report  presents  an  overview  of  a  body  of  theoretical  work  that  has  been  the  foundation  for 
several  laser  guide-star  experiments  conducted  by  Lincoln  Laboratory  [3-5].  Although  the  original  intent 
of  this  investigation  was  to  develop  a  means  to  propagate  near-diffraction-limited  laser  beams  through  the 
atmosphere,  recent  changes  in  classification  guidelines  have  made  it  possible  to  apply  these  techniques 
to  ground-based  astronomical  telescopes.  In  keeping  with  the  original  emphasis  on  near-perfect  turbu¬ 
lence  compensation  for  short  wavelengths  and  marginal  seeing  conditions,  however,  the  actuator  densities 
and  correction  bandwidths  discussed  are  generally  much  higher  than  those  recommended  by  other  authors 
[6]. 

To  lay  the  proper  groundwork  for  the  design  guidelines  presented  at  the  end  of  this  report,  it  is 
necessary  to  first  review  the  inherent  limitations  of  the  phase-conjugation  process  and  justify  the  need  for 
synthetic-beacon  turbulence  probes.  Section  2  of  this  report  concentrates  on  conventional  adaptive-optics 
systems;  in  this  context,  the  term  “conventional”  refers  to  systems  having  access  to  a  real  beacon  located 
within  the  isoplanatic  patch  of  the  target  object.  Section  2  also  introduces  the  turbulence-  and  wind-profile 
models  that  are  the  basis  for  all  the  numerical  calculations  that  follow.  The  isoplanatic  errors  associated 
with  the  deployment  of  synthetic  beacons  are  the  main  topic  of  Section  3,  where  a  brief  comparison  of 
single-  and  multiple-beacon  geometries  is  also  presented. 

The  main  thrust  of  this  investigation  is  presented  in  Section  4,  in  which  a  methodology  is  proposed 
for  balancing  the  major  sources  of  compensation  error  to  optimize  the  performance  of  an  adaptive-optics 
system.  Recognizing  that  constructs  exploiting  natural-star  beacon  illumination  are  highly  desirable  from 
a  cost  perspective,  Section  4  contains  a  quantitative  discussion  of  fractional  sky  coverage  for  both  natural- 
and  synthetic-beacon  systems.  In  conclusion,  significant  benefits  accrue  with  the  use  of  Irser  guide  stars, 
but  the  development  of  more  compact  and  less  expensive  laser  sources  will  be  necessary  before  this 
approach  becomes  practical  for  widespread  application. 
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2.  FUNDAMENTAL  CONCEPTS 


To  avoid  ambiguity  in  the  systems  analyses  that  follow,  a  brief  discussion  of  some  of  the  more 
important  turbulence  issues  is  given  in  this  section.  Of  particular  interest  is  the  treatment  of  Strehl  ratio 
(on-axis  inter.,  ty  relative  to  a  diffraction-limited  image)  and  resolution  for  systems  capable  of  partial 
phase  compensation. 


2.1  STATISTICAL  DESCRIPTION  OF  ATMOSPHERIC  TURBULENCE 

A  quantitative  treatment  of  optical  turbulence  usually  begins  with  a  statistical  description  of  the 
index  of  refraction.  A  useful  quantity  is  the  structure  function,  which  represents  the  expected  variance  in 
refractive  index  between  two  points 

A.s(("(n)-«(^))2}  ,  (2-1) 

where  and  r2  are  three-dimensional  vectors.  Kolmogorov  [7]  has  shown  that  in  the  inertial  subrange 
this  function  is  isotropic  and  proportional  to  the  2/3  power  of  the  scalar  distance,  r, 

0„(r)  =  C„2«)rw  .  (2-2) 

The  Cl  parameter  characterizes  the  strength  of  the  refractive  index  fluctuations  as  a  function  of  the 
altitude,  h.  From  this  expression.  Fried  [8]  derived  the  form  of  the  phase  structure  function  at  the  surface 
of  a  collection  aperture 

Z)0(r)  =  2.9U2r5/3sec(OjC„2(/t)rf/i  ,  (2-3) 

where  k  =  2n/k  is  the  optical  wavenumber,  f  is  the  zenith  angle,  and  the  integral  is  performed  over  the 
beam  path.  (MKS  units  have  been  assumed  in  computing  the  leading  constant.)  Fried  found  it  convenient 
to  define  a  turbulence  coherence  length,  r0,  over  which  the  gross  wavefront  distortion  can  be  described 
as  a  uniform  tilt.  Thus, 


^(r)-6.88(r/rQ)5/3  ,  (2-4) 

where 

ic 

rQ  =  {o.423*2  sec(£)  J  C2(/i)<//r}~  *  (2'5) 

An  aperture  of  dimension  r0  will  produce  a  nearly  diffraction-limited  far-field  image  that  will  appear  to 
change  position  as  the  atmosphere  evolves  with  time.  This  parameter  establishes  the  critical  spatial 
^’mension  of  a  subaperture,  within  which  a  unique  phase  error  can  be  measured  and  subsequently 
corrected  by  the  movement  of  a  single  mechanical  actuator. 
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2.2  TURBULENCE  AND  WIND- VELOCITY  MOMENTS 


Turbulence  calculations  make  widespread  use  of  turbulence-profile  integrals  of  the  type  indicated 
in  Equation  (2-3).  The  n01  moment  is  defined  in  the  following  manner: 

K.  =/£?(<)<■<«  ,  (2-6) 

where  the  integral  is  taken  over  the  entire  turbulence  region.  The  low-order  moments  are  most  sensitise 
to  ground-level  turbulence,  whereas  the  high-order  moments  are  more  heavily  influenced  by  upper- 
altitude  turbulence.  Using  this  notation,  the  expression  for  the  coherence  length  would  be 

r0  =  {0.423/t2sec(O^0}"3/5  *  (2'5>) 

The  isoplanatic  patch,  which  defines  the  angle  within  which  the  distortion  over  the  turbulence  path  will 
be  essentially  unchanged,  can  be  shown  to  depend  on  the  5/3  moment 

«0  =  {2.9U2sec8^O|<^(^)A5/3^}~3/5={2.9U2sec8/3(O//53}~35  '  (2'7) 

The  temporal  characteristics  of  turbulence-induced  phase  perturbations  are  often  modeled  on  the 
assumption  that  turbulence  can  be  described  as  a  set  of  frozen  layers,  each  of  which  moves  across  the 
aperture  at  speeds  that  vary  as  a  function  of  altitude.  For  simplicity,  the  direction  of  motion  is  assumed 
to  be  perpendicular  to  the  plane  containing  the  zenith  angle.  The  velocity  moments  have  the  following 
form: 


*’"(<)<«  ,  (2-8) 

where  v(h )  is  the  wind-velocity  function.  A  critical  time  constant  specifying  the  interval  over  which 
turbulence  remains  essentially  unchanged  derives  from  Greenwood  [9],  and  can  be  expressed  as  a  func¬ 
tion  of  the  5/3  velocity  moment 

tq  =|2.91^2sec(^)Jc2{A)r'5^3(A)aK|  =  j2.91£2sec(£’)j'5^3j  •  ^-9) 

To  attach  numerical  values  to  these  critical  system  parameters,  baseline  profiles  must  be  established 
for  turbulence  and  wind  velocity.  During  the  past  two  decades,  several  dozen  models  have  been  devel¬ 
oped  in  an  attempt  to  describe  refractive-index  measurements  made  at  a  variety  of  locations  throughout 
the  world.  For  daytime  conditions  at  inland  sites,  the  SLC-Day  (standardized  under  the  Submarine-Laser 
Communications  program)  and  the  Hufnagel-Valley  turbulence  models  are  frequently  used  for  systems 
analysis.  [See  Equations  (C-l)  and  (C-2)  in  Appendix  C.}  To  the  best  of  our  knowledge,  however,  there 
are  no  standard  turbulence  profiles  for  good  astronomical  sites  and  nighttime  viewing  conditions  where 
the  seeing,  defined  as 
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seeing 


(2-10) 


4 


is  often  better  than  0.5  arc-seconds  [10].  Therefore,  the  authors  have  elected  to  assemble  a  modified 
version  of  the  Hufnagel-Valley  profile  that  yields  a  20-cm  r0  and  a  20-prad  0O  for  0.55  pm  radiation  at 
zenith 


C;(ft)  =  8.16xl0"54  /t10  exp(-/i/1000)  +  3.02xl0~17  exp(-A/1500) 

(2- II) 

+  1.90xl0-15  exp(~/i/100) 

Note  that  the  altitude  is  given  in  meters  and  the  units  of  Cn2  are  m'm. 

Temporal  fluctuations  in  the  cumulative  phase  distortion  are  somewhat  less  critical  to  this  analysis; 
therefore,  for  the  wind  profile  a  standard  form  proposed  by  Bufton  is  used; 

v(k)  =  5  +  30  exp|-[(/i  -  9400)/4800]2  }  (2.l2) 

The  units  for  this  expression  are  m/sec. 

The  turbulence  and  wind  models  are  plotted  as  a  function  of  altitude  in  Figures  1  and  2,  respec¬ 
tively.  A  comparison  of  key  turbulence  parameters  can  be  found  in  Table  1. 
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Figure  1.  Comparison  of  the  HV-21  model  typically  used  for  calculations  of  daytime  turbulence  and  the  modified 
HV  model  developed  by  the  authors  to  represent  nighttime  turbulence  conditions  at  good  seeing  locations. 
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Figure  2.  The  Bufton  wind-velocity  profile. 
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TABLE  1 

Summary  of  Turbulence  Parameters  for  Imaging  Systems  Operating  at  0.55  pm 


ro 

(cm) 

8»e«ing 

(arc*sec) 

So 

(mrad) 

*0 

(msec) 

SLC-Day 

Model 

0°  Zenith 

D 

n 

13.8 

n 

45°  Zenith 

D 

3.0 

8.0 

2.2 

HV-21 

Model 

0°  Zenith 

5.6 

2.5 

7.8 

2.1 

45°  Zenith 

Q 

3.1 

a 

1.7 

Modified  HV 
Model 

0°  Zenith 

20.0 

0.69 

20.0 

45°  Zenith 

16.2 

0.85 

11.5 

5.1 

23  PERFORMANCE  MODELS  FOR  CONVENTIONAL  ADAPTIVE-OPTICS  SYSTEMS 

An  adaptive-optics  system  for  real-time  phase  conjugation  can  be  characterized  as  a  highly-parallel 
servo  device  capable  of  sensing  and  correcting  the  incoming  wavefront  at  the  pupil  plane  of  an  imaging 
sensor.  Performance  close  to  the  diffraction  limit  of  the  input  aperture  can  be  achieved  in  the  limit  that 
the  angular  separation  between  the  turbulence  probe  and  the  target  object  is  small  (<0()),  the  spacing 
between  control  elements  on  the  active  optical  element  is  well  matched  to  the  turbulence  coherence  length 
(<r0),  and  a  sufficiently  high  update  rate  is  maintained  (>l/r0).  In  most  practical  implementations  of  this 
concept,  however,  performance  compromises  are  necessary  to  reduce  component  costs  and  maximize  the 
sensitivity  of  the  wavefront  sensor.  At  this  point,  a  brief  description  of  the  architecture  of  an  adaptive- 
optics  system  is  appropriate.  A  more  detailed  discussion  of  the  basic  components  and  an  interesting 
overview  of  the  early  history  of  this  field  can  be  found  in  John  Hardy’s  classic  paper  on  this  subject  [11]. 

The  four  essential  elements  of  a  conventional  adaptive-optics  system  are  shown  in  Figure  3.  Light 
originating  from  an  exoatmospheric  source  is  corrupted  in  both  amplitude  and  phase  as  a  result  of  random 
fluctuations  in  the  refractive  index  of  the  intervening  air.  The  quality  of  the  image  formed  at  the  focus 
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of  a  ground-based  telescope  is  largely  driven  by  phase  distortions,  which  can  be  corrected  through  the 
insertion  of  an  optical  surface  having  the  conjugate  optical  path  difference  (OPD).  Most  compensation 
systems  apply  this  conjugate  phase  with  a  pair  of  active  elements  consisting  of  a  high-speed  tilt  mirror 
and  a  deformable  mirror  that  removes  the  figure  (tilt-removed)  component  of  the  distortion.  Error  signals 
are  generated  by  the  phase  sensor,  which  actually  measures  the  first  derivative  of  the  phase  of  the 
incoming  wavefront.  A  phase-computation  device,  referred  to  as  a  reconstructor,  transforms  the  output 
of  the  wavefront  sensor  into  a  set  of  drive  signals  that  control  the  active  optical  elements. 

Each  of  these  components  has  a  unique  set  of  erTor  mechanisms  that  together  establish  the  overall 
effectiveness  of  the  adaptive-optics  system.  In  the  following  sections,  the  first-order  effects  of  these  errors 
are  quantified;  the  most  important  effects  are  then  combined  in  a  servo-loop  model  that  provides  a  concise 
description  of  the  dynamic  performance  characteristics. 


Figure  3.  Essential  components  of  an  adaptive -optics  phase-compensation  system.  Error  signals  developed  from 
the  phase  difference  between  the  incoming  wavefront  and  the  deformable  mirror  are  measured  by  the  phase  sensor 
and  subsequently  applied  to  the  two  active  optical  elements. 


23.1  The  Wavefront  Sensor 

Essentially  all  wavefront-sensing  instruments  in  current  use  measure  local  changes  in  the  pupil- 
plane  phase  function;  a  direct  determination  of  phase  is  difficult  because  a  coherent  full-aperture  reference 
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would  be  required.  At  least  three  sensing  techniques  have  been  successfully  demonstrated;  Hartmann 
sensors,  in  which  the  input  pupil  is  segmented  by  a  set  of  lenslets  that  form  separate  images  of  the  source 
in  the  Fourier  plane  [11];  shearing  sensors  that  convert  phase  gradients  into  intensity  variations  [11];  and 
sensors  that  measure  local  wavefront  curvature  [12].  The  noise  characteristics  of  the  Hartmann  and 
shearing  sensors  are  similar  and  reasonably  well  understood,  although  the  former  is  more  easily  de¬ 
scribed.  For  this  reason,  the  Hartmann  design  will  be  used  as  the  prototype  in  this  analysis. 

The  sensor  model  used  in  this  study  is  based  on  the  realization  illustrated  in  Figure  4.  A  lenslet 
array  placed  at  an  image  of  the  pupil  plane  is  used  to  generate  a  set  of  independent  far-field  images;  the 
positions  of  these  images  are  linearly  related  to  the  wavefront  tilt  over  the  lenslets'  diameter.  An  array 
of  sensors  is  then  employed  to  measure  the  centroid  position  of  each  of  the  focused  spots. 

Although  the  first  detector  arrays  were  composed  of  discrete  photomultiplier  tubes,  these  were  later 
replaced  by  image  intensifies  mated  to  solid-state  cameras  and,  more  recently,  by  high-quantum-effi- 
ciency  CCDs  [13].  Because  the  noise  level  of  a  CCD  camera  is  proportional  to  the  square  root  of  the  read¬ 
out  bandwidth,  it  is  often  advantageous  to  use  separate  imagers  to  perform  the  x-  and  y-gradient 
measurements.  This  configuration  allows  on-chip  data  summation  to  be  performed  in  the  direction  per¬ 
pendicular  to  the  gradient  measurement  and  reduces  the  number  of  data  points  from  n2  to  2 n  per  subaperture. 
Under  these  conditions,  a  one-dimensional  sensor  model  incorporating  a  linear  subarray  of  n  pixels  is 
sufficient  to  describe  all  relevant  effects.  The  size  of  the  subarray  will  depend  on  the  required  dynamic 
range,  although  a  pair  of  elements  per  subaperture  per  camera  should  be  adequate  for  continuous  wave 
(CW)  operation  under  benign  turbulence  conditions. 

To  simplify  this  discussion  of  phase-sensor  performance,  it  will  be  assumed  that  the  lenslets  divide 
the  pupil  plane  into  an  array  of  subapertures  of  dimension  d}.  (The  physical  components  are  generally 
much  smaller  than  the  collection  aperture  and  would  be  placed  in  a  collimated  region  of  the  optical  path 
at  the  output  of  the  telescope.)  For  a  phase  change  of  one  wave  across  the  subaperture,  the  Hartmann  spot 
will  experience  a  linear  translation  of  A  f/ds,  where  /  is  the  lenslet  focal  length  and  A  is  the  wavelength 
of  the  turbulence  probe.  Near  the  null  position  of  a  two-element  centroid  detector,  an  appropriate  linear 
algorithm  for  determining  the  local  gradient  would  be 


g  =  2it 


-1 


-rL  +  rR 


2A/r(0)J  r^  +  r* 


(rad  /  subaperture) 


(2-13) 


9 


20731** 


Figure  4.  Schematic  description  of  a  Hartmann  wavefront  sensor.  The  incoming  wavefront  is  divided  into  a  matrix 
of  subapertures  by  a  lenslet  array,  which  produces  a  set  of  focused  spots  onto  a  detector  array.  The  displacement 
of  the  spots  from  their  local  null  position  provides  a  measure  of  the  local  phase  gradient. 


where  II  x)  is  the  normalized  intensity  distribution  of  the  Hartmann  spot  and  the  symbols  TL  and  rR 
represent  the  number  of  photoelectrons  collected  by  the  left  and  right  detector  elements,  respectively.  If 
the  spot  is  diffraction-limited. 


r(x) 


ds  >jfsin(g</Jx/A/)|2 

A/J  1  Xdsx/kf  J 


(2-14) 


for  a  square  lenslet,  where  I TO)  is  simply  the  inverse  of  the  spot  diameter,  dspoV  When  the  subaperture 
becomes  much  larger  than  r0,  the  spot  diameter  will  be  roughly  proportional  to  Xf/rQ.  Yura  [14]  suggests 
the  following  profile  for  the  general  case: 
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Combining  Equations  (2-13)  and  (2-15)  yields  the  following  result: 


g  =  1ty  1  + 


Y 2 


1 


Kr0  J 


~TL  +  rR 
rL  +  rR 


(2-16) 


The  photon  noise  associated  with  the  centroid  detector  characterized  by  Equation  (2-16)  is  readily 
developed  by  applying  Poisson  statistics  to  the  detector  outputs,  so  that  the  count  variance  in  each  pixel 
is  equal  to  the  absolute  number  of  photoelectrons.  Using  N =  2 (/}_  +  7^)  to  represent  the  total  number 
of  photoelectrons  collected  by  the  x  and  y  sensors  within  a  single  subaperture,  an  expression  for  the 
gradient  variance  for  a  beam  near  the  null  position  is  derived: 


(2-17) 


The  second  major  error  is  a  result  of  random  signal  fluctuations  generated  by  the  transfer  of 
photoelectrons  through  the  detector  array,  off-chip  preamplifiers,  and  digitizers.  This  uncorrelated  noise 
variance,  represented  by  N2^,  can  be  reduced  relative  to  the  photon-noise  contribution  if  an  optical 
intensifier  having  a  net  photoelectron  gain  of  Ge  is  placed  at  the  entrance  aperture  of  the  wavefront 
sensor.  Assuming  that  each  pixel  has  an  identical  noise  contribution,  the  following  expression  is  obtained 
from  Equation  (2-16): 


(2-18) 


A  more  general  expression  for  subapertures  that  incorporate  more  than  two  elements  is  derived  in  Section 
C.2  of  Appendix  C. 

Equations  (2-17)  and  (2-18)  provide  an  accurate  description  of  the  performance  of  a  Hartmann 
sensor;  to  maintain  generality,  however,  a  few  additional  concepts  will  be  introduced  at  this  time.  To 
allow  for  the  possibility  that  the  compensation  wavelength,  Ac,  differs  substantially  from  the  phase  sensor 
wavelength,  Ap,  the  multiplicative  factor  [kc/kp^  must  be  inserted  to  convert  the  error  expression  to  rad2 
of  phase  error  at  Xc.  Furthermore,  because  it  is  convenient  to  define  the  fundamental  turbulence  param¬ 
eters  (rQ,  0O,  and  T0)  at  the  compensation  wavelength,  it  can  be  seen  from  Equation  (2-5)  that  the 
turbulence  coherence  length  at  Xp  is  properly  given  by  (k//kc)6l5r0.  Finally,  the  following  representation 
will  be  applied  for  the  number  of  collected  photoelectrons: 


N  = 

pe 


2  n 


hck 


n  <£  4 


p  j 


(2-19) 
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where  lp  (W/m2)  is  the  irradiance  at  the  entrance  of  the  wavefront  sensor,  r?  is  the  sensor  quantum 
efficiency,  rd  is  the  sensor  dwell  time,  and  d2  is  the  subaperture  collection  area.  The  final  result  is 
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(2-20) 

Notice  that  the  photon-noise  term  is  independent  of  the  subaperture  size  when  ds  is  much  larger  than  the 
coherence  length  for  Xp.  In  this  regime,  the  noise  reduction  expected  with  an  increased  subaperture 
diameter  (due  to  the  larger  number  of  detected  photons)  is  directly  offset  by  a  reduction  in  sensitivity  to 
wavefront  gradients. 

The  motivation  for  using  low-noise  CCD  arrays  as  opposed  to  intensified  cameras  for  wavefront 
sensing  is  addressed  in  Figure  5.  Although  the  gain  of  an  intensifier  can  always  be  made  high  enough 
to  obtain  photon-noise-limited  operation,  the  quantum  efficiency  of  these  devices  seldom  exceeds  10 
percent.  In  contrast,  backside-illuminated  CCD  cameras  are  now  available  that  are  better  than  80  percent 
efficient  in  the  visible  and  can  operate  at  high  frame  rates  with  transfer  noise  levels  in  the  20-40-electron 
range.  Within  the  normal  range  of  operation  for  high-quality  turbulence  compensation  (where  the  average 
gradient  error  would  be  lower  than  A/10),  the  performance  of  a  well-designed  CCD  system  is  expected 
to  surpass  that  of  the  intensified  camera.  For  this  reason,  the  wavefront  sensor  in  current  use  at  Lincoln 
Laboratory  incorporates  a  pair  of  custom-built  64  X  64-CCD  focal  planes. 


2C7319-5 


Figure  5.  Performance  comparison  of  the  two  major  types  of  Hartmann-sensor  cameras.  This  calculation  was  per¬ 
formed  for  a  phase  sensor  having  two  pixels  per  subaperture  subarray  and  a  subaperture  dimension  smaller  than  ra 
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2 32  The  Phase  Reconstructor 

In  the  limit  of  infinite  measurement  resolution,  the  relationship  between  the  phase,  <p,  and  the  phase 
gradient,  g ,  is  given  by  the  Poisson  equation,  with  Neumann  boundary  conditions, 

V20  =  div  g  .  (2-21) 

For  a  discrete-measurement  system,  this  relationship  can  be  represented  by  a  linear  process 

g  =  A0  ,  (2-22) 

in  which  g  and  0  are  vectors  and  A  is  a  two-dimensional  matrix.  The  least-square  solution  to  this 
equation  can  be  recovered  from  the  generalized  inversion  of  the  A  matrix 

.  (2-23) 

Because  B  is  invariant  for  a  given  sensor  configuration,  the  matrix  inversion  only  needs  to  be  performed 
once  and  can  be  preprogrammed  into  the  computational  device. 

Modem  reconstructors  employ  digital  processors  that  can  perform  the  matrix  operation  indicated 
in  Equation  (2-23)  within  a  time  period  much  smaller  than  the  turbulence  time  constant.  For  example, 
the  reconstructor  built  by  Lincoln  Laboratory  for  a  set  of  recent  field  experiments  performs  a  241  X  436 
matrix  calculation  with  8-bit  precision  in  less  than  50  psec;  this  corresponds  to  an  average  computational 
rate  of  the  order  of  109  operations  per  second  [15].  Mathematical  techniques  are  currently  under  devel¬ 
opment  that  would  permit  much  larger  arrays  to  be  processed  even  more  efficiently. 

The  noise  characteristics  of  a  reconstructor  are  defined  by  its  error  propagator,  which  is  a  multi¬ 
plicative  constant  that  relates  the  variance  of  the  input  gradient  measurements  to  the  variance  of  the 
output  phase  estimates.  Herrmann  [16]  has  shown  that  the  error  propagator  for  a  compensation  system 
that  has  a  few  hundred  actuators  is  approximately  unity. 

In  summary,  a  well-designed  reconstruction  system  will  introduce  no  significant  temporal  delay  or 
excess  noise  into  the  phase-compensation  process. 

233  The  Deformable  Mirror 

Several  innovative  mirror  designs  have  been  developed  in  the  last  ten  years  that  incorporate  either 
modal  (Zemike  deconvolution)  [12]  or  zonal  control  of  the  mirror  surface  [17,18];  both  continuous- 
surface  and  segmented-mirror  designs  have  been  successfully  demonstrated.  The  14-cm  deformable  mirror 
in  current  use  at  Lincoln  Laboratory  is  an  Itek  design  that  incorporates  a  set  of  241  discrete  actuators 
mated  to  a  thin  glass  facesheet,  as  illustrated  in  Figure  6.  The  characteristics  of  the  facesheet  have  been 
chosen  to  provide  a  smooth  transition  between  actuator  drive  points  under  normal  turbulence  conditions 
and  to  yield  a  fitting  error  of  better  than  1/20  wave  rms  over  a  60-cm  telescope  aperture.  The  actuators 
are  made  by  stacking  thin  sections  of  lead  magnesium  niobate  (PMN),  which  is  a  low-voltage  electrostrictive 
material.  A  deflection  of  ±1.5  pm  has  been  measured  with  a  drive  voltage  of  ±150  volts;  the  average  full- 
stroke  settling  time  is  on  the  order  of  200  psec.  At  this  speed,  the  temporal  characteristics  of  the  mirror 
have  a  negligible  effect  on  the  overall  system  performance. 
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Figure  6.  Cross-section  of  a  typical  deformable  m  '.  • ror  employing  a  flexible  facesheel  and  a  two-dimensional  array 
of  actuator  stacks.  The  actuators  are  cemented  to  the  facesheet  so  that  the  surface  can  be  deformed  by  either 
extending  or  contracting  the  actuator  dimension. 


The  major  degrading  effect  associated  with  a  deformable  mirror  is  fitting  error,  which  occurs 
whenever  the  actuator  spacing  in  the  pupil  plane  exceeds  the  turbulence  coherence  length.  Greenwood 
[  1 9]  has  derived  a  relationship  between  the  error  variance  and  the  number  of  actuators,  Na  ; 


K-)«„g-a274(D/r»)5,3N;5/6 


which  can  be  rewritten  to  reflect  the  error  associated  with  a  subaperture  spacing  of  ds  : 


(2-24) 


(2-24’) 


The  general  form  of  this  equation  has  been  verified  in  simulation  studies  using  the  measured  parameters 
of  Lincoln’s  deformable  mirror,  although  the  leading  constant  has  been  found  to  vary  between  0.4  and 
0.8,  depending  on  how  the  outer  edge  of  the  flexible  mirror  surface  is  constrained.  These  effects  can  be 
minimized  by  employing  at  least  one  ring  of  active  actuators  outside  the  illuminated  region  of  the  aperture 
to  provide  a  smooth  transition  between  the  mirror  interior  and  the  support  collar.  For  the  purpose  of  this 
study,  the  following,  somewhat  more  conservative,  estimate  will  be  adopted: 


figure 


(2-25) 


23.4  A  Time-Dependent  Model  of  the  Phase-Compensation  Process 

The  individual  characterizations  of  component  performance  provide  no  information  about  the  dy¬ 
namic  effects  that  drive  the  update  rate  of  an  adaptive-optics  system.  This  issue  is  perhaps  best  addressed 
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by  constructing  a  servo-loop  model  for  each  of  the  active  control  elements.  Thus,  two  classes  of  models 
will  be  required — one  to  describe  the  movement  of  the  actuators  on  the  deformable  mirror  and  one  for 
the  control  of  the  tilt  mirror.  The  final  result  will  be  a  set  of  error  expressions,  the  first  of  which  describes 
the  total  figure  distortion: 


K“"Lv„=  +  jWF/,(f)df  (rad2  of  phase) ,  (2-26) 

.0  .0 

'  V  '  V  -v  . '' 

turbulence -induced  figure  error  figure  noise 

and  similarly  for  the  overall  tilt  error, 

oo  oo 

K-L.'  iWWfW  +  (rad2  of  tilt).  (2-27) 

0  0 

v - - - -  - - ' 

turbulence -induced  tilt  error  tilt  noise 

The  first  term  in  each  of  these  equations  is  the  residual  error  due  to  uncompensated  turbulence,  whereas 
the  second  term  is  the  noise  associated  with  a  finite  signal  into  the  sensor.  The  parameters  T(f)  and  F(f) 
icpresent  servo  transmission  functions  and  signal  power  spectra,  respectively.  Because  the  power  spectra 
for  figure  and  tilt  are  quite  different,  the  bandwidth  requirements  for  each  must  be  developed  separately. 

A  schematic  diagram  of  the  servo  system  is  provided  in  Figure  7.  The  wavefront  sensor  measures 
the  difference  between  the  optical  phase  of  the  incoming  wavefront  and  the  phase  applied  to  the  deformable 
mirror.  This  error  signal  is  averaged  over  a  dwell  time,  Xj,  and  subsequently  passed  to  an  accumulator 
that  adds  the  measured  error  to  the  previous  mirror-drive  position.  For  simplicity,  the  mirror  update  is 
assumed  to  occur  as  a  step  input  at  the  end  of  each  integration  cycle,  and  each  mirror  correction  is  held 
until  the  next  update.  This  quasi-CW  representation  accurately  describes  the  operation  of  modem  digital- 
control  architectures. 
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Figure  7.  Servo-loop  diagram  for  a  digital  control  system  employing  a  boxcar  averaging  device  and  a  digital 
accumulator.  The  error  signal  at  the  wavefront  sensor  is  averaged  for  a  period  of  Xd  The  accumulator  provides 
the  deformable  mirror  with  an  impulse  drive  signal  a’  the  end  of  each  integration  interval;  this  signal  is  subse¬ 
quently  held  until  the  next  drive  signal. 


Frequency-domain  models  for  each  of  the  components  shown  in  Figure  7  can  be  found  in  a  number  of 
standard  texts  [20].  The  phase  error  sampling  system  is  described  by  the  function 


»,(/)* 


sin(nfrd) 


x  exp(— iVr/rrf)  = 


l-exp(— i2itfrd) 
2nfrd 


(2-28) 


which  includes  an  integration  period  of  xd  and  an  average  delay  time  of  r d/2  between  the  input  and  output 
signals.  The  digital  accumulator,  which  adds  the  incremental  error  signal  to  the  current  mirror  position, 
is  represented  by 


«.(/)* 


8i 

2rcfTd  * 


(2-29) 
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where  St  is  a  normalized  loop-gain  constant.  Finally,  the  drive  response  for  a  mirror  time  constant,  xm, 
would  be 


illtf  Tm  +1  ‘ 


(2-30) 


Generally,  x _«t.  for  state-of-the-art  hardware,  so  this  term  can  be  safely  ignored. 


From  the  component  transfer  functions  given  above,  an  equivalent  control-circuit  description  can 
be  developed  for  the  propagation  of  optical  turbulence  and  electrical  noise  through  the  system.  The  net 
loop  response  to  atmospheric  turbulence  is  found  to  be 


WH 


(2nftdy 


\  +  Hs(f)Ha(f)\ 


4si n2(nfrd) 


(2  nfrdf 


-gt  +(2  jcfrdy 


(2-31) 


The  noise  transfer  function  obtains  from  the  assumption  that  electrical  noise  originates  in  the  sensing 
device  and  is  transferred  to  the  control  elements  through  the  signal  averager  and  accumulator.  Therefore, 


rm  I  „ _ ! _ 

+  «,(/)»„(/)  ,  (Ixftrf  t  (2 K/T,)4 

gt  4s?s  m2(nfrd) 


(2-32) 


These  two  equations  are  plotted  as  a  function  of  normalized  frequency  and  gain  in  Figures  8  and  9.  The 
transfer  curve  for  turbulence  has  the  characteristics  of  a  high-pass  Filter,  whose  shape  most  closely 
matches  an  ideal  step  function  when  the  gain  is  approximately  equal  to  one.  By  contrast,  the  noise  transfer 
curve  for  this  system  is  a  low-pass  filter,  which  is  also  optimized  when  g(  =  1 .  Both  of  these  expressions 
are  greatly  simplified  for  gain  of  unity,  resulting  in  the  approximations 


T  (/)  (27r/T^)2  +  (2?r/Trf)4/12 

*  l  +  {2nfrd)A/l2 


(2-33) 


Uf)- 


l  +  (2xfTd)*/l2 


(2-34) 


The  turbulence-rejection  bandwidth  of  the  first  of  these  functions  is  approximately  1/10 r^,  whereas  the 
integral  of  the  noise  transfer  function  is  about  1/3t^ 


z 

o 


o 

z 

3 

u. 

a 

UJ 

u. 

(A 

Z 

< 

oc 

t- 


NORMALIZED  FREQUENCY  (2?rtrd) 


Figure  8.  Turbulence  transfer  function  for  a  quasi-CW  closed-loop  compensation  system  as  a  function  of  loop  gain. 
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Figure  9.  Noise  transfer  function  for  a  quasi-CW  closed-loop  compensation  system  as  a  function  of  loop  gain. 


18 


To  complete  this  analysis,  four  sets  of  signal  spectra  are  required  to  describe  the  effects  of  turbu¬ 
lence  and  noise  on  both  the  figure-  and  tilt-compensation  systems.  The  figure  component  of  turbulence 
is  developed  from  a  paper  by  Greenwood  and  Fried  [21],  which  is  discussed  in  more  detail  in  Section 
C.3  of  Appendix  C.  This  spectrum,  which  is  plotted  in  Figure  10,  can  be  approximated  by  a  function 
having  two  segments 


'0.132sec(C)^2D4^J2/5r"73,5/4/3  :/<  0.705  D*1^35^5 

0.0326sec(0*c  *>5/3  /~8'3  -f  *  0.705 *  (2-35) 

The  units  are  rad2  of  phase  error  per  Hz.  A  concise  treatment  of  full-aperture  tilt  was  recently  published 
by  Tyler  [22].  This  spectrum  also  has  a  simple  form  in  its  low-  and  high-frequency  limits  that  can  be 
joined  at  their  intersection  to  give 


1.60  sec(f)  t»_w  fm  :/<  0.445  D"'  «C$  *, |4’3 

0.028  scc{()D-svmr'w  :/>0.U5D-‘v:jg^so  (2-36) 


where  the  units  here  are  rad2  of  single-axis  tilt  per  Hz  and  normalized  for  integration  over  positive  values 
of  /.  A  plot  of  this  function  for  a  4-m  aperture  is  given  in  Figure  1 1 . 


Figure  JO.  Power  spectrum  for  turbulence-induced  figure  error.  This  model  incorporates  atmospheric  moments 
developed  from  the  modified  Hufhagel-Valley  turbulence  profile  and  the  Bufton  wind  profile. 
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Figure  II.  Power  spectrum  for  turbulence-induced  tilt  error,  given  in  units  of  rad2  of  single-axis  tilt  per  Hz  This 
model  incorporates  atmospheric  moments  developed  from  the  modified  Hufnagel-  Valley  turbulence  profile  and  the 
Bufton  wind  profile. 


The  figure  component  of  the  noise  power  spectrum  can  be  obtained  from  Equation  (2-20)  by 
assuming  that  the  noise  is  white  and  passes  through  a  filter  having  the  form 

H  (/)  = 

avgVJ  f 

which  is  the  transfer  function  for  a  boxcar  averager  with  a  time  constant,  rd.  The  integral  of  this  transfer 
function  introduces  a  normalization  factor  of  lx#  so  that 

F/h(/)  =  2Td(<Tfigure)noise 


sin(*/TJ 


(2-37) 
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A  Hartmann  sensor  can  also  be  viewed  as  a  tracking  device  in  which  a  phase  shift  of  1  rad  over  a  receiver 
diameter  D  corresponds  to  a  1/kD  angular  displacement  of  the  focal  spot.  Thus  the  noise  spectrum  for 
tracking  can  be  obtained  from  Equation  (2-38)  by  converting  the  collection  area  to  -f  D2,  the  wavenumber 
to  kt,  the  dwell  time  to  Tdt,  the  irradiance  to  /,,  and  multiplying  by  (1/fc.D)2; 
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(2-39) 


An  evaluation  of  all  the  individual  error  sources  is  now  possible  by  combining  the  last  six  equations 
and  performing  the  integrations  indicated  in  Equations  (2-26)  and  (2-27).  The  noise  terms  are  straight¬ 
forward,  because  the  integrals  of  the  figure  and  tilt  transfer  functions  are  1/3^  and  l/3r^,  respectively. 
The  turbulence  integrals  are  somewhat  more  complicated,  but  are  amenable  to  a  few  simplifying  assump¬ 
tions.  For  high-bandwidth  figure  compensation,  accurate  results  are  obtained  by  using  the  high-frequency 
form  of  the  figure  spectrum  over  the  entire  frequency  range,  so  that 

]T,(f)F,(f)df  - 0.0326sec(fKJ J  (2k/z<)  f-mdf 

o  0  l  +  (2/r/  id)  /12 

=  2.80sec(£)h2  v5p  3  (rad2  of  phase)  .  (2-40) 

Therefore,  the  figure  error  associated  with  the  operation  of  the  servo  system  is  found  to  be 
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The  total  figure  error  is  obtained  by  adding  the  deformable  mirror  fitting  error  [specified  in  Equa¬ 
tion  (2-25)]  to  this  result. 
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The  compulation  of  the  turbulence-induced  tilt  error  is  performed  by  first  recognizing  that  the  high- 
frequency  roll-off  in  the  tilt  spectrum  is  sufficiently  rapid  to  allow  contributions  beyond  the  low-fre¬ 
quency  cut-off  to  be  ignored.  In  this  range,  the  transfer  function  is  proportional  to  the  square  of  the 
frequency.  From  Equations  (2-33)  and  (2-36), 


0  445  D~'  vils 

ao  v  r_(Jj  Vj4/3 

J  T,(f)F,(f)if  -1.60  sec(f )  „_l/3  f  ( 2<lf  td,  f  rmdf 

0  0 

«  4. 09sec(£) zr7/3  ^  t^Jf  T2dt  ( rad2  of  tilt) ,  (2-42) 

which  leads  to  the  following  expression  for  the  total  single-axis  tilt  jitter: 
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Equations  (2-41)  and  (2-43)  provide  a  reasonably  complete  and  accurate  description  of  the  dynamic 
characteristics  of  an  adaptive-optics  servo  system. 


2.4  WORKING  DEFINITIONS  OF  ERROR  VARIANCE,  STREHL,  AND  IMAGING 
RESOLUTION 


Quantitative  estimates  of  overall  system  performance  are  notoriously  difficult  to  achieve,  particu¬ 
larly  when  the  cumulative  errors  are  large.  The  usual  approach  is  to  compute  the  phase  error  in  the  pupil 
plane  due  to  each  source,  then  develop  a  total  error  by  adding  the  variances 


(2-44) 


This  will  result  in  an  over-estimate  of  the  error  unless  the  individual  effects  are  statistically  independent. 


For  laser-beam  propagation,  the  quantity  of  interest  is  the  Strehl  ratio,  which  is  defined  as  the 
relationship  between  the  measured  on-axis  intensity  and  the  peak  intensity  of  a  diffraction-limited  beam. 
It  is  observed  that  when  the  rms  pupil-plane  phase  distortion  is  much  smaller  than  one  wave,  the  width 
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of  the  short-exposure  beam  will  remain  essentially  unchanged  but  the  peak  energy  will  degrade  exponen¬ 
tially  according  to  the  expression 


Strehl  *  exp^-tJ2) 


(2-45) 


which  is  known  as  the  extended  Marshal  approximation.  The  energy  lost  from  the  main  beam  will 
appear  as  a  set  of  diffraction-limited  “speckles”  that  are  randomly  distributed  over  a  region  that  is 
inversely  proportional  to  the  turbulence  coherence  diameter.  As  the  distortion  level  is  increased,  the  main 
beam  will  eventually  become  so  depleted  that  it  becomes  indistinguishable  from  the  rest  of  the  speckle 
field.  In  this  regime  it  is  more  difficult  to  establish  an  unambiguous  definition  of  Strehl,  but  the  form 
suggested  by  Yura  [14]  is  probably  appropriate; 


Strehl  - - j 

l  +  [o/P0] 


(2-46) 


p0  is  a  characteristic  coherence  length  that  can  be  computed  for  either  short  or  long  exposures;  for 
uncorTected  turbulence  both  of  these  quantities  are  approximately  equal  to  rQ. 

No  single  strategy  for  specifying  Strehl  or  sensor  resolution  is  suitable  for  all  situations.  In  many 
cases,  however,  a  reasonably  accurate  estimate  of  these  quantities  can  be  obtained  by  combining  the 
Strehl  expressions  given  by  Mardchal  and  Yura,  under  the  assumption  that  the  energy  lost  from  the  main 
peak  is  distributed  over  a  beam  profile  characterized  by  A JpQ.  The  resulting  beam  profile  for  short 
exposure  would  be 


[Strehl  ]SE  =  exp(-<7f,gUre)+ 
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l  +  [D/pff 


(2-47) 


2 

where  the  short-exposure  coherence  length,  ps^,  may  also  depend  on  CTf,gUre.  From  this  description  of 
the  Strehl,  a  form  for  the  composite  system  resolution  is  evolved  that  incorporates  a  weighted  average 
of  the  diffraction-limited  resolution,  1.22A C/D,  and  the  resolution  in  the  limit  of  large  phase  error. 


(l.22AcH/l  +  [D/pff. 


The  following  expression  is  the  rms  of  the  product  of  the  Strehl  and  the 


resolution  for  the  two  limiting  cases: 


[Resolution]SE 


(2-48) 


An  equivalent  equation  for  the  long-exposure  Strehl  can  be  developed  by  increasing  the  short- 
exposure  beam  profile  to  reflect  the  addition  of  a  single-axis  tilt  jitter  (rads  of  angle)  of  .  Recognizing 


23 


that  a  Gaussian  beam  characterized  by  a  standard  deviation  of  0.45 X/D  accurately  reproduces  the  core 
of  the  Airy  function,  it  can  be  seen  that 
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and 


[Resolution]LE 


(2-50) 


where  the  tilt  error  for  the  large-phase  error  limit  has  been  included  in  the  long-exposure  coherence 
length,  p gE.  According  to  Yura  [14],  the  long-exposure  coherence  length  is  approximately  equal  to  r0; 


(2-51) 


and  the  short-exposure  coherence  length  is  always  somewhat  larger; 


1  +  0-37  (r0/O)’/3 

»  - 


(2-52) 


The  application  of  this  somewhat  unorthodox  approach  to  beam  characterization  can  be  demon¬ 
strated  in  a  short  discussion  of  fitting  error  effects  resulting  from  low-resolution  phase  conjugation. 
Several  investigators  have  claimed  that  near  diffraction-limited  resolution  can  be  achieved  with  only  a 
few  actuators  [6],  arguing  that  the  long-exposure  beam  profile  will  always  consist  of  a  sharp  central  peak 
surrounded  by  a  diffuse  background.  This  statement  is  valid  as  long  as  the  tracking  system  can  reliably 
discriminate  between  the  main  lobe  of  the  beam  and  the  competing  speckles  in  the  short-exposure  image. 
Consider,  for  example,  the  two  short-exposure  profiles  shown  in  Figure  12,  which  represent  the  output 
of  a  ray-trace  simulation  incorporating  a  compensation  system  in  which  the  ratio  of  the  actuator  spacing 
to  the  turbulence  coherence  length  was  changed  from  1 .8  to  2.5.  In  the  latter  case,  the  transition  out  of 
the  Marshal  regime  has  begun;  it  is  much  more  difficult  to  unambiguously  locate  and  track  the  central 
lobe.  As  indicated  in  Figure  13,  the  effective  short-exposure  beamwidth  is  clearly  turbulence-dominated 
as  soon  as  the  phase  error  exceeds  4  rad2  (ds/rQ*=  3.5).  Note  that  the  simulation  results  are  in  good 
agreement  with  Equation  (2-48). 
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Figure  12.  Comparison  of  short-exposure  beam  profiles  near  the  transition  to  strong  phase  distortion.  This 
simulation  was  performed  for  an  aperture  diameter  of  4  m  and  an  ro  of  16  cm.  The  upper  figure  represents  a 
d/ro  of  1.8,  for  which  a  Strehl  of  0.2  was  obtained.  The  lower  figure  corresponds  to  a  d/r0  of  2.5  and  a  Strehl  of 
0.03.  (The  vertical  scale  for  the  lower  figure  has  been  multiplied  by  a  factor  of  5.) 
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Figure  13.  Comparison  of  the  predicted  short-exposure  resolution  with  far-field  beam-profile  measurements  de¬ 
rived  from  a  ray-trace  simulation;  Kolmogorov  phase  screens  were  compensated  with  a  zonal-type  adaptive  element 
having  a  variable  actuator  spacing.  Each  simulation  result  represents  the  average  over  5  independent  screens  that 
were  each  adjusted  to  produce  an  rQ  of  16  cm  over  a  4-m  aperture. 
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3.  LASER  GUIDE-STAR  CONCEPTS 


The  efficacy  of  adaptive  optics  has  been  repeatedly  demonstrated  in  laboratory  and  field  experi¬ 
ments  conducted  within  the  defense  community  over  the  last  15  years.  At  this  point  there  is  little  doubt 
that  well-compensated  images  can  be  obtained  with  ground-based  telescopes,  provided  that  the  object  is 
bright  enough  to  supply  a  few  hundred  photoelectrons  to  each  subaperture  per  dwell  time.  For  a  nominal 
set  of  conditions  (a  16-cm  subaperture,  5-msec  dwell  time,  and  a  20  percent  throughput),  this  corresponds 
to  an  input  flux  of  order  107  photons/m2-sec,  which  is  roughly  equivalent  to  a  magnitude  7  object  for 
visible  observations.  This  limitation  would  obviously  exclude  many  objects  of  interest  to  astronomers. 

The  concept  of  using  synthetic  sources  for  turbulence  sensing  was  first  proposed  to  the  high- 
energy-laser  beam-control  community  in  1982  [1].  Shortly  thereafter,  a  classified  program  was  initiated 
to  investigate  the  theory  and  implementation  of  adaptive-optics  devices  incorporating  artificial  beacons. 
Lincoln  Laboratory  was  one  of  several  organizations  involved  in  that  effort  and,  to  the  best  of  our 
knowledge,  obtained  the  first  experimental  evidence  of  laser  guide-star  compensation  in  the  summer  of 
1988  [4].  By  the  spring  of  1991,  Strehl  ratios  exceeding  0.3  were  routinely  obtained  in  the  visible  with 
a  60-cm  telescope  and  a  241 -actuator  adaptive-optics  system.  Figure  14  shows  an  example  of  pre-  and 
post-compensation  data  of  the  star  Procyon.  In  this  experiment  the  wavefront  sensor  was  driven  by  a 
beacon  produced  by  a  512-nm  dye  laser,  the  return  from  which  was  range-gated  at  6  km. 

Within  the  astronomy  community,  an  independent  research  effort  has  been  active  since  1985  [2], 
and  in  the  last  few  years  the  literature  on  this  subject  has  expanded  considerably  [23-31].  Although  many 
of  the  important  system  design  considerations  have  already  been  discussed  in  these  articles,  the  previous 
treatments  have  been  largely  qualitative.  This  section  will  provide  a  quantitative  treatment  of  the  perfor¬ 
mance  limitations  that  arise  from  a  variety  of  anisoplanatic  effects. 

3,1  FOCAL  ANISOPLANATIC  ERRORS  PRODUCED  BY  A  SINGLE  ARTIFICIAL  BEACON 

Anisoplanatic  errors  will  occur  in  any  sampling  geometry  in  which  the  measured  turbulence  path 
does  not  coincide  with  the  path  followed  by  the  image  radiation.  As  indicated  in  Figure  15,  light  rays 
originating  from  a  distant  source  propagate  along  parallel  lines  to  the  receiver  aperture.  The  measurement 
of  accumulated  phase  distortion  afforded  by  a  laser  guide  star  is  inherently  imperfect  as  a  result  of 
unsampled  turbulence  above  the  beacon  and  incorrectly  sampled  turbulence  below  the  beacon.  This  form 
of  error  is  referred  to  as  focal  anisoplanatism,  because  it  results  from  a  difference  in  range  between  the 
beacon  and  the  celestial  object. 
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Figure  14.  Short-exposure  visible-wavelength  images  of  Procyon  taken  in  February  iyw  wim  the  60-cm  telescope 
operated  by  Lincoln  Laboratory  at  Maui.  The  upper  picture  is  a  1-msec  exposure  taken  with  a  CCD  camera  having 
a  19-prad  field  of  view.  The  lower  image  was  collected  a  millisecond  later  after  the  phase  distortion  had  been 
corrected  with  a  24 1 -actuator  deformable  mirror.  A  512-nm  synthetic  probe  placed  at  an  altitude  of  6  km  was  used 
to  measure  the  turbulence-induced  phase  distortion.  The  estimated  Strehl  ratios  for  the  uncompensated  and  com¬ 
pensated  images  are  0.05  and  0.4,  respectively. 
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Figure  15.  Illustration  of  the  two  major  sources  of  phase  error  introduced  by  laser  guide-star  sources.  The  vertical 
rays  represent  light  originating  from  a  distant  source  that  accumulates  phase  error  in  traveling  through  the 
atmosphere.  The  beacon  rays  follow  slightly  different  paths  to  the  receiver  and  are  unable  to  sample  the  turbulence 
above  the  beacon  or  to  properly  sample  the  turbulence  below  the  beacon. 


An  analytic  description  of  phase  errors  associated  with  anisoplanatic  effects  is  most  easily  obtained 
by  developing  the  appropriate  ensemble  averages  of  the  phase  structure  function  in  the  pupil  plane  of  the 
receiver.  An  approach  developed  by  Sasiela  [32]  icads  to  error-variance  expressions  having  the  following 
general  form: 


o2  =  0.207 k2c  Jcn2(z){J/(K-)g(^z)dic}d2  f  (3-D 

0 

where  the  z  integration  is  performed  over  the  slant  range  from  the  receiver  to  the  object  and  the  k 
integration  is  taken  over  the  entire  two-dimensional  space  of  the  pupil  plane,  /(v)  is  the  turbulence 
spectrum,  which,  for  most  problems,  is  assumed  to  be  the  normalized  Kolmogorov  spectrum 

f(K)=  iCun  .  (3-2) 
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The  turbulence-sensing  geometry  is  characterized  by  the  filter  function.  g(ic,c),  and  will  be  unique  to 
each  problem.  As  an  example,  consider  the  residual  phase  error  for  a  system  in  which  the  aperture- 
averaged  tilt  is  removed.  The  filter  function  that  eliminates  the  average  phase  (piston)  and  the  x-  and  v- 
tilt  over  an  aperture  of  diameter  D  is 

\2J,(KD2)f  r4y,(vD2)f 


‘  [  kD‘  2  J  [  i 

from  which  the  following  expression  is  obtained: 

0  0  <- 
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The  figure  error  associated  w'ith  the  uncorrected  turbulence  above  the  beacon  is  derived  directly 
from  Equation  (3-4)  by  setting  the  lower  limit  of  the  integral  over  z  to  the  beacon  altitude,  H.  If  the  n**1 
upper  turbulence  moment  is  defined  to  be 

oo 

,  (3-5) 

H 

the  error  due  to  the  unsampled  turbulence  can  be  written  as 

<4*,  =°-057  D5,3*t2sec(CK(W)  (3.6) 

The  phase  error  below  the  beacon  arises  from  the  difference  between  the  parallel  rays  traveling  between 
the  scattering  altitude  and  the  ground  and  the  rays  originating  from  the  beacon.  In  this  case,  the  filter 
function  is  somewhat  more  complicated  and  leads  to  the  expression 


kD:2 


kD(H~z)I2H 
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The  solution  is  most  easily  expressed  as  a  series  expansion  involving  the  low-altitude  turbulence  mo¬ 
ments,  computed  from  ground  level  to  the  beacon.  An  accurate  representation  can  be  obtained 

with  the  first  two  expansion  terms 


rr2  ~  D513 

°lowcr  ~  U 


K :  sec(£) 
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rim 
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rim 
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(3-8) 


Note  that  both  the  low-  and  high-altitude  equations  exhibit  an  explicit  D'/3  k2  sec(£)  dependence  that 
is  decoupled  from  the  beacon  altitude  and  turbulence  profile  dependencies  contained  in  the  expansion 
series. 


The  predicted  focal  anisoplanatic  error  for  visible  imaging  with  a  4-m  aperture  is  plotted  as  a 
function  of  beacon  altitude  in  Figure  16.  A  single  beacon  placed  at  10  km  should  reduce  the  figure  error 
to  about  3  rad2,  which  (based  on  the  analysis  of  resolution  shown  earlier  in  Figure  13)  is  just  within  the 
range  in  which  a  unique  diffraction-limited  central  lobe  would  be  formed.  A  90-km  source  would  yield 
an  image  having  a  much  higher  signal-to-clutter  ratio  and  would  produce  a  Strehl  ratio  of  approximately  0.6 


Figure  16.  Comparison  of  phase  variance  contributions  due  to  turbulence  lying  above  and  below  the  synthetic 
beacon.  The  shape  of  the  low-altitude  curve  is  complicated  by  the  inclusion  of  a  larger  fraction  of  the  atmosphere 
as  the  beacon  is  raised. 
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22  IMPLEMENTATION  AND  PERFORMANCE  OF  MULTIPLE-BEACON  SYSTEMS 


The  efficacy  of  systems  employing  a  single  synthetic  beacon  came  into  question  fairly  early  in 
Lincoln’s  investigation  of  guide-star  techniques  and  subsequently  led  to  the  development  of  a  variety  of 
alternate  beacon  geometries.  Because  the  figure  error  due  to  focal  anisoplanatism  has  been  shown  to  have 
a  D5/3  dependence,  there  is  a  potential  benefit  in  deploying  multiple  beacons,  each  of  which  would  be 
devoted  to  correcting  one  section  of  the  complete  aperture.  This  concept  is  illustrated  in  Figure  17.  If  a 
•otal  of  Nb  artificial  sources  were  to  be  employed,  the  effective  collection  aperture  for  the  residual 
turbulence  below  the  beacons  would  be  reduced  by  a  factor  of  jV^12,  so  that 


flower  =ty> 5/6  £5/3*c2sec(£) 
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Unfortunately,  the  practical  implementation  of  this  scheme  is  far  from  trivial,  and  the  following 
issues  are  noted. 


1)  In  order  to  mate  the  aperture  sections  to  their  respective  beacons,  a  discrimination  approach  using 
time  delay,  color,  polarization,  or  fteld-of-view  must  be  employed  so  that  the  signals  from  the  wavefront 
sensor  can  be  correctly  processed.  If  the  beacons  are  launched  sequentially  in  time,  the  entire  pattern  must 
be  generated  within  the  atmospheric  time  constant. 

2)  Because  wavefront  sensors  measure  phase  gradients  rather  than  absolute  phase,  a  means  must  be 
developed  to  measure  gradients  across  the  seams  between  the  aperture  sections.  To  accomplish  this 
“stitching”  process,  the  detectors  at  the  edges  of  the  section  seams  must  be  capable  of  measuring  the  light 
from  both  of  the  neighboring  sources.  If  the  beacons  are  launched  in  p?ra!lel,  this  may  necessitate  the 
use  of  multiple  sensors  with  overlapping  ftelds-of-view. 
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Figure  17.  Multiple-beacon  sampling  geometry.  Focal  anisoplanatism  is  reduced  by  positioning  a  beacon  over  the 
center  of  each  mirror  section.  "Stitching’'  errors  occur  when  the  relative  source  positions  cannot  be  accurately 
determined  and  result  in  low  spatial-frequency  figure  distortions. 
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The  section-stitching  issue  has  been  studied  at  some  length  within  the  beam-control  community,  and  there 
is  still  considerable  doubt  as  to  the  ability  of  multiple-beacon  arrays  to  produce  results  that  are  signifi¬ 
cantly  better  than  those  achievable  with  single-beacon  systems.  The  basic  problem  arises  from  the  need 
to  accurately  fix  the  relative  positions  of  the  sources  at  the  scattering  layer.  Figure  17  illustrates  the  ideal 
geometry  in  which  each  of  the  beacons  lies  directly  above  the  center  of  its  respective  section.  If  one  of 
the  spots  were  to  be  laterally  displaced  within  this  pattern,  however,  all  the  subapertures  within  its  section 
would  measure  an  anomalous  tilt  that  would  result  in  a  low  spatial-frequency  distortion  of  the  recon¬ 
structed  phase.  Such  displacements  can  be  expected  to  occur  as  a  result  of  imperfections  in  the  beacon- 
projection  system;  on  a  more  fundamental  level,  fluctuations  in  position  are  unavoidable  due  to  the  the 
fact  that  each  beam  is  projected  through  a  slightly  different  turbulence  path.  The  use  of  overlapping 
gradient  measurements  can  reduce  this  effect,  but  there  will  always  be  a  non-correctable  error  component. 

The  uncorrectable  part  of  the  stitching  error  is  best  described  in  the  context  of  a  measurement 
process  that  attempts  to  map  the  actual  beacon  positions  by  allowing  each  section  to  simultaneously  view 
every  source.  In  the  absence  of  turbulence,  this  mapping  process  would  be  error-free.  When  the  atmo¬ 
sphere  is  present,  however,  non-zero  tilt  measurements  resulting  from  source  displacements  become 
confused  with  turbulence-induced  tilt  differences. 

A  detailed  evaluation  of  the  multiple-beacon  phase-reconstruction  process  shows  that  the  net  re¬ 
sidual  is  related  to  the  difference  between  the  section  tilts  obtained  from  the  artificial  beacon  array  and 
the  target  object.  To  perform  this  calculation,  the  section-tilt  variance  due  to  focal  anisoplanatism  must 
be  determined  and  the  cumulative  effect  of  this  error  over  an  aperture  having  a  large  number  of  sections 
must  be  estimated.  The  first  of  these  calculations  is  performed  by  applying  the  filter  function  for  aperture- 
tilt  to  Equation  (3-1).  For  a  section  having  a  diameter,  £>,  the  appropriate  filter  function  is 


tilt  phase 


4J2(kDs/2) 

kDs/2 


(3-10) 


which  leads  to  the  following  expansion  for  the  tilt  error  due  to  turbulence  below  the  beacon: 


4,Ph«-0f  *c2sec({)|o.368i^+-J 

An  evaluation  of  the  net  error  over  the  full  aperture  requires  knowledge  of  the  correlation  statistics  of 
these  fluctuations.  (If  the  section-tilt  errors  were  perfectly  correlated,  no  high  spatial-frequency  distortions 
would  be  produced  and  the  figure  error  would  be  unaffected.)  This  correlation  matrix  has  been  evaluated, 
and  in  Section  B.3  of  Appendix  B  the  error  propagator  for  a  system  having  a  large  number  of  beacons 
is  shown  to  be  approximately  equal  to  0.1 1(D/Z))5/3.  Therefore,  the  stitching  component  of  the  figure 
error  is 


Etching  =05/3*f2sec(O  0.040 


Aw 


(3-12) 
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which  has  the  same  Z>5/3  k 3  sec(£)  dependence  found  in  the  previous  error-variance  expressions  and  is 
independent  of  the  number  of  sections. 

As  suggested  by  Equation  (3-9),  it  is  possible  in  principle  to  reduce  the  low-altitude  anisoplanatism 
term  to  a  negligible  level  by  deploying  a  sufficient  number  of  synthetic  beacons.  In  the  limiting  case,  the 
major  contributors  to  the  total  error  would  be  upper-altitude  anisoplanatism  and  section-stitching  noise. 
This  result  is  shown  in  Figure  18,  and  in  comparison  with  Figure  16  it  can  be  seen  that  a  significant  error 
reduction  is  predicted  for  all  beacon  altitudes.  However,  this  improvement  must  be  carefully  balanced 
against  the  increased  complexity  of  the  laser  system,  wavefront  sensor,  and  phase  reconstructor. 
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Figure  18.  Figure  error  expected  for  the  combination  of  unsampled  upper-uliitude  turbulence  and  aperture-section 
stitching  with  multiple  beacons. 


The  practical  implementation  of  the  stitching  concepts  just  described  has  been  demonstrated  by 
Lincoln  Laboratory  in  a  four-beacon  laboratory  simulation  performed  in  1988  and  in  a  two-beacon  field 
experiment  in  November  1990  [5].  In  the  latter  experiment,  a  pair  of  pulsed  dye  lasers  was  used  to 
generate  sources  that  were  separated  in  time  by  approximately  1  msec  and  displaced  laterally  by  half  the 
diameter  of  the  aperture.  In  both  of  these  tests,  a  small  improvement  in  the  far-field  Strehl  was  recorded 
when  the  multiple-oeacon  images  were  compared  with  the  single-beacon  results 
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33  LONG-EXPOSURE  TRACKING  CONSIDERATIONS 


All  synthetic-beacon  concepts  suffer  from  a  common  deficiency — because  the  laser  radiation  fol¬ 
lows  the  same  turbulence  path  in  both  the  transmit  and  backscatter  directions,  no  information  can  be 
gained  about  the  probe  position  relative  to  the  target  object.  Unless  the  error  due  to  turbulence-induced 
tilt  jitter  can  be  sensed  and  eliminated,  the  long-exposure  image  obtained  with  guide-star  figure  correction 
will  be  no  better  than  that  produced  by  an  uncompensated  receiver. 

For  imaging  engagements  in  which  exposures  much  longer  than  the  atmospheric  time  constant  are 
required,  it  will  be  necessary  to  stabilize  the  imaging  system  on  an  exoatmospheric  fiducial  object  that 
is  within  the  isoplanatic  angle  for  overall  tilt.  Recall  that  the  standard  definition  for  isoplanatic  angle  is 

d0={29lkc  sec8/3(0^5/3}  ’  (3'I3) 


which  is  based  on  the  following  description  of  the  phase  error  resulting  from  two  beam  paths  separated 
by  an  angle  &. 


beacon  offset 


(3-14) 


Sasiela  [32]  has  derived  an  expression  for  the  error  variance  associated  with  overall  tilt,  where  in  this  case 
the  aperture  filter  function  is 
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This  expression  represents  the  product  of  three  terms:  the  first  is  a  conversion  from  phase  variance  to 
single-axis  angular  tilt,  the  second  is  the  tilt  filter  function  introduced  in  Equation  (3-10),  and  the  third 
is  the  function  for  angular  offset.  The  resulting  error  (in  units  of  rad2  of  single-axis  tilt)  is  approximately 
given  by 
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where  Hc  =  DjO  sec(£)  is  the  altitude  at  which  the  two  beam  paths  no  longer  overlap.  The  first  term 
of  this  series  is  usually  adequate  as  long  as  Hc  is  greater  than  10  km. 

Combining  the  first  term  of  Equation  (3-16)  with  Equation  (2-50)  yields  the  long-exposure  reso¬ 
lution  for  tilt  anisoplanatism: 


[Resolution  ]LE  =  1.22 


(l) 

i 

M 

2 

,  1  +  4.94 

{ 

1 

A 

V  c ) 

\ 

(<&)■ 


beacon  offset 


where  the  isoplanatic  angle  for  tilt  is  defined  to  be 


=  {0.668^ sec3(^2D-^3}  1/2  .  (3-18) 

For  a  set  of  nominal  system  parameters  (D  =  4  m,  A.  =  0.55  pm,  and  £  =  45°)  and  the  modified  HV 
turbulence  model,  (fig)^  =  40  prad;  this  number  is  only  a  factor  of  3.5  larger  than  60  for  these  conditions. 

A  plot  of  the  long-exposure  resolution,  developed  from  Equations  (3-16)  and  (3-17),  is  given  in 
Figure  19.  This  curve  displays  a  linear  dependence  on  the  separation  angle  when  is  between  1 

and  6,  but  for  larger  values  of  6  an  asymptotic  limit  is  reached  that  is  V2  larger  than  the  single-beam 
tilt  jitter  for  uncompensated  turbulence. 

This  research  indicates  that  a  guide-star  compensation  system  must  comprise  two  separate,  and 
equally  essential,  servo  elements — a  deformable  mirror  controlled  by  the  laser  beacon  and  a  tracking 
mirror  slaved  to  a  bright  fiducial  star.  The  probability  of  finding  a  suitable  fiducial  object  within  the  tilt 
isoplanatic  angle  will  be  addressed  in  the  next  section. 


Figure  19.  Resolution  as  a  function  of  angular  displacement  between  the  target  object  and  the  fiducial  star.  This 
function  is  approximately  linear  between  50  and  250  prod. 


36 


4.  OPTIMIZATION  AND  PERFORMANCE  ANALYSIS  FOR  A  4-m  SYSTEM 


At  this  point,  an  optimized  telescope  design  may  be  proposed  and  its  performance  characteristics 
analyzed.  To  motivate  the  need  for  synthetic  beacons,  brightness  requirements  will  also  be  developed  for 
celestial  objects. 

To  simplify  this  analysis,  a  baseline  system  configuration  will  be  used  for  all  calculations  that 
follow.  The  list  of  parameters,  given  in  Table  2,  includes  a  4-m  aperture,  a  viewing  angle  of  45°,  and 
a  common  wavelength  of  0.55  pm  for  wavefront  sensing,  tracking,  and  imaging.  The  sensor  character¬ 
istics  of  an  80  percent  quantum  efficiency  and  20  noise  electrons  are  nominal  for  high-ffame-rate  CCD 
camera  systems. 

Note  also  that  a  total  throughput  parameter,  y,  will  appear  in  all  of  the  error  expressions  in  this 
section  so  that  the  computed  irradiance  values  can  be  referenced  to  the  top  of  the  atmosphere.  Estimates 
of  throughput  can  vary  greatly  and  will  depend  on  the  complexity  of  the  collection  system  and  the  manner 
in  which  the  light  is  apportioned  between  the  various  sensors.  For  the  purpose  of  this  study  a  value  of 
20  percent  will  be  used  for  both  the  phase  sensor  and  tracking  systems. 

Unless  otherwise  stated,  the  parameters  listed  in  Table  2  should  be  assumed  for  all  numerical 
results. 

4.1  OPTIMIZATION  OF  SUBAPERTURE  DIMENSION  AND  DWELL  TIME 

The  three  most  important  sources  of  error  for  conventional  (non-guide-star)  adaptive-optics  systems 
have  been  shown  to  be  fitting  error,  phase-sensor  noise,  and  finite-bandwidth  error.  By  combining  Equa¬ 
tions  (2-25)  and  (2-41),  the  residual  figure  distortion  for  a  system  that  senses  wavefront  phase  at  a 
wavelength  Xp  and  images  at  wavelength  Ac  is  obtained.  Thus, 
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TABLE  2 

Baseline  System  Parameters  for  Optimization  Study 


Aperture  Diameter 

D  *  4  m 

Wavelength 

Ap  =  A,  =  Ac  =  0.55  pm 

Zenith  Angle 

;  =  45° 

Tcia:  System  Throughput 

=  Y,  =  0.2 

Detector  Quantum  Efficiency 

.3 

It 

o 

bo 

Detector  Noise  Level 

"rms  =  20 

r0  -  16.2  cm 

60  =  11.5  grad 

Turbulence  Parameters 

r0  =  5.1  msec 

(0o)tm  *  40-3  Mfad 
(Xo)^  =  12.3  msec 

Pq  *  2.65  x  10~13m1/3 

Turbulence  Moments 

=  1.78  x  10~7mz 

P2  =  3.96  x  10-®m7/3 

u_1/3  =  1.40  x  10r13sec1/3 

Wind  Velocity  Moments 

=  1.23  x  10"11m2-sec5/3 
u14/3  =  2.28  x  10~7m5-sec14/3 

The  optimization  of  the  subaperture  dimension  and  phase  sensor  dwell  time  is  simplified  if  the  third  term 
in  Equation  (4-1)  is  rewritten  for  the  limit  ds»rQ.  A  general  treatment  of  this  problem  is  given  in  Section 
C.4  of  Appendix  C,  and  from  Equation  (C-20)  it  can  be  shown  that  the  relative  contributions  of  the  three 
error  terms  appearing  in  Equation  (4-1)  are  y^-,  y^,  and  yf,  respectively.  For  a  specified  figure-variance 
design  goal,  the  following  set  of  system  characteristics  is  developed: 
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(4-4) 


The  last  of  these  equations  can  be  rewritten  to  show  the  explicit  dependence  on  the  number  of  phutcelec 
trons  collected  by  each  subaperture 


(4-4') 

These  expressions  indicate  that  ds/r0  and  Td/x0  will  typically  be  close  to  unity  and  that  the  number  of 
photoelectrons  that  must  be  collected  by  each  subaperture  per  measurement  interval  will  be  on  the  order 
of  a  few  hundred.  A  plot  of  ds  and  xd  as  a  function  of  the  allowed  Figure  error  is  given  in  Figure  20. 
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Figure  20.  Optimization  analysis  of  the  phase-compensation  system  parameters  as  a  function  of  figure  variance. 
A  total  throughput  value  of  20  percent  has  been  assumed  in  this  calculation. 
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The  tracking-sysfem  dwell  time  can  be  optimized  in  a  similar  fashion,  although  in  this  case  it  will 
be  assumed  that  the  bandwidth  of  the  tilt  sensor  is  low  enough  to  permit  photon-noise-limited  operation. 
From  Equation  (2-43), 
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(rad2  of  single-axis  tilt),  (4-5) 


where  a  critical  time  constant  for  tilt  has  been  defined  following  the  procedure  used  earlier  to  develop 
the  tilt  isoplanatic  angle: 


(4-6) 


When  the  dwell  time  is  properly  adjusted,  ‘he  relative  contributions  of  the  two  error  terms  in  Equation 
(4-5)  will  be  -j  and  y,  respectively.  Thus, 
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represents  the  optimal  setting,  where  the  jitter  radius  has  been  normalized  to  the  half-width  of  the 
diffraction-limited  beam  at  the  compensation  wavelength.  The  tracker  irradiance  requirement  is 
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and  at  this  level  the  detector  will  collect  the  following  number  of  photoelectrons  within  a  single  measure¬ 
ment  cycle: 
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Once  again  the  dwell  time  must  be  roughly  matched  to  the  atmospheric  time  constant  if  good  performance 
is  to  be  achieved  and  the  requisite  photoelectron  count  is  of  the  order  of  a  few  hundred. 
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4.2  PERFORMANCE  EXPRESSIONS  FOR  THE  OPTIMIZED  SYSTEM  DESIGN 

By  reformulating  the  optimization  equations  to  express  them  in  terms  of  source  irradiance  and 
combining  them  with  their  respective  error  expressions,  a  compact  set  of  variance  relationships  is  ob¬ 
tained; 
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Quantitative  estimates  of  the  long-exposure  Strehl  and  resolution  follow  from  Equations  (2-49)  and 
(2-50),  which  can  be  simplified  if  the  system  is  assumed  to  operate  in  the  range  in  which  the  extended 
Marechal  approximation  is  valid  (°figurc<  3  rad2),  with  the  result 
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[Resolution]LE  =  1.22^^-j  ^1  +  4.94 


(2-50’) 


A  plot  of  these  functions  for  lp-  lt  is  given  in  Figure  21.  The  equivalent  stellar  magnitude  is  also  given 
for  wideband  sensing  in  the  visible,  where  the  relationship 

mv  =  -2.5  log(/)- 21.2  (4-11) 


has  been  applied. 


Figure  21.  Performance  analysis  of  the  optimized  figure -  and  tilt-compensation  systems  as  a  function  of  irradiance 
(referenced  to  the  top  of  the  atmosphere).  A  total  throughput  value  of  20  percent  has  been  assumed  for  both  systems. 


43  SKY-COVERAGE  CALCULATIONS  FOR  NATURAL  SOURCES 

The  conclusion  that  can  be  drawn  from  Figure  21  is  that  phase  compensation  at  visible  wavelengths 
requires  relatively  bright  objects  (mv<  9),  whereas  high-precision  tracking  can  be  performed  with  stars 
much  dimmer  than  16th  magnitude.  In  this  context  the  advantage  of  guide-star  compensation  becomes 
obvious,  as  a  natural  source  is  only  required  to  close  the  tracking  loop.  In  order  to  quantify  the  expected 
improvement,  it  will  be  useful  to  obtain  estimates  of  the  fraction  of  the  sky  that  can  be  covered  with 
natural  and  synthetic  beacons.  The  fractional  sky  coverage  is  defined  to  be 

Fractional  Sky  Coverage  =  it  t?2  x  Average  Star  Density  ,  (4- 1 2) 

where  t?  represents  a  radius  that  can  be  drawn  around  each  bright  star  within  which  the  error  due  to  tilt 
anisoplanatism  will  be  tolerable  according  to  some  set  of  system  performance  goals. 

Star  density  computations  have  been  performed  based  on  information  extracted  from  the  Infrared 
Handbook  [33],  which  gives  average  population  estimates  as  well  as  numbers  for  the  galactic  pole  and 
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equator.  The  resulting  density  curves  as  a  function  of  visual  magnitude  are  given  in  Figure  22.  For 
computational  purpose  it  will  be  convenient  to  replace  these  curves  with  analytical  functions;  the  fol¬ 
lowing  expressions  provide  an  excellent  fit  for  visual  magnitudes  less  than  20: 


Average  Density  =  1.45  exp(0.96  m^) 

(stars/rad2)  , 

(4-13) 

Density  at  Pole  =  1.27  X  lO^m^  2 

(stars/rad2)  , 

(4-14) 

Density  at  Equator  =  3.97  exp(mv) 

(stars/rad2)  . 

(4-15) 

Figure  22.  Comparison  of  stellar  densities  for  the  galactic  pole,  the  galactic  equator,  and  the  average  density  as 
a  function  of  visual  magnitude. 


Equations  (4-11)  and  (4-13)  can  be  combined  to  produce  a  first-order  estimate  of  the  average  density  as 
a  function  of  inradiance: 

Average  Density  =  2.02xl0-9  7"104  (stars/rad2) .  (4-16) 

The  area  of  coverage  afforded  by  each  bright  star  depends  on  the  level  of  error  due  to  tilt 
anisoplanatism  that  can  be  tolerated.  From  Equations  (3-13)  and  (3-14)  it  can  be  seen  that  to  maintain 
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a  specified  figure  error  of  (^figure) 
than 


beacon  offset 


,  the  angular  displacement  from  the  beacon  must  be  less 
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This  equation  is  applicable  when  natural  stars  are  used  as  the  light  source  for  both  figure  and  tilt 
compensation.  When  a  synthetic  beacon  is  used  to  drive  the  wavefront  sensor,  the  fiducial  star  must  be 


close  enough  to  the  target  object  to  sustain  a  tilt  jitter  of  less  than 
(3-16)  and  (3-18), 
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4.4  LASER  REQUIREMENTS  FOR  GUIDE-STAR  COMPENSATION 


The  final  component  of  this  discussion  is  an  estimate  of  the  power  required  to  generate  a  laser  guide 
star.  The  expressions  for  atmospheric  backscatter  are  well  known,  and  good  tutorials  on  this  subject  can 
be  found  in  texts  by  Hinkley  [34]  and  Measures  [35]. 

In  Figure  16  strong  evidence  is  presented  for  the  use  of  sodium  beacons  to  achieve  high-quality 
phase  compensation  at  visible  wavelengths.  To  fully  exploit  resonance  emission,  however,  the  laser  line 
width  must  be  well  matched  to  the  D2  absorption  line  (-3  GHz),  and  saturation  effects  resulting  from  the 
long  natural  decay  time  of  the  excited  state  (-16  nsec)  must  be  avoided.  The  second  of  these  parameters 
results  in  a  limit  of  about  30  kW/m2  on  the  peak-power  density  of  the  projected  beam  [30]. 

Saturation  can  be  reduced  by  either  increasing  the  length  of  the  laser  pulse  or  the  size  of  the  laser 
spot  within  the  scattering  region.  In  Lincoln’s  theoretical  studies,  virtually  no  degradation  of  the  quality 
of  the  phase  measurement  was  found  as  long  as  the  projected  beam  was  smaller  than  the  aperture 
diameter  (or  the  section  diameter  for  a  multiple-beacon  system).  To  avoid  any  degradation  of  the  phase 
sensor’s  sensitivity,  however,  the  beacon  must  not  be  so  large  that  it  becomes  resolved  by  a  subaperture. 
For  most  scenarios  of  interest,  this  latter  restriction  will  establish  the  practical  upper  limit  on  the  laser- 
spot  diameter. 

Assuming  a  laser  pulse  width  of  200  psec  and  a  50-cm  beacon  diameter,  the  limit  on  the  pulse 
energy  would  be  approximately  1  J.  This  number  is  slightly  greater  than  the  energy  requirement  that  is 
derived  below  for  CW  operation  in  the  visible;  therefore,  saturation  effects  can  be  neglected  for  most 
first-order  calculations. 

For  quasi-CW  operation  using  a  high  rep-rate  beacon  laser,  the  backscattered  irradiance  from 
sodium  layer  measured  at  the  entrance  to  the  wavefront  sensor  is  predicted  to  be 
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where 


Ip  -  average  irradiance  at  the  wavefront  sensor  (W/m2), 

=  average  laser  power  (W), 

7 1  —  throughput  of  the  laser  beam  to  the  scattering  layer  =  0.2, 

£  =  zenith  angle  =  45°, 

Hs  -  sodium  layer  altitude  =  9  X  104  m, 

Cs  =  sodium  column  abundance  =  5  x  1013  atoms/m2,  and 
da{  n)/dQ  =  scattering  cross-section  =  6.6  x  10~17  m2/sr. 

The  combination  of  Equations  (4-4)  and  (4-19)  produces  an  expression  for  the  average  nower  required 
for  CW  phase  compensation; 


(4-20) 

A  plot  of  the  average  laser  power  and  pulse  repetition  rate  as  a  function  of  the  allowable  figure  variance 
is  given  in  Figure  23.  A  total  two-way  throughput  of  4  percent  was  assumed  in  the  development  of  these 
laser  requirements. 
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Figure  23.  Beacon-laser  specifications  for  resonance  backscatter  front  the  earth’s  sodium  layer.  A  two-way 
throughput  of  4  percent  has  been  assumed  in  this  calculation. 


The  development  of  sodium-frequency  lasers  has  been  active  at  Lincoln  Laboratory  for  several 
years.  The  most  recent  efforts  have  focused  on  a  solid-state  design  based  on  the  sum-frequency  mixing 
of  two  diode-pumped  Nd:YAG  lasers  operating  at  1.064  pm  and  1.319  pm  [36].  A  24-mJ  840-Hz  laser 
was  constructed  in  1990,  and  it  is  predicted  that  this  same  technology  could  be  used  to  build  a  device 
capable  of  generating  several  hundred  watts  of  average  power. 

4.5  PERFORMANCE  COMPARISONS  FOR  NATURAL-STAR  ANT*  LASER  GUIDE-STAR 
SYSTEMS 

All  of  the  relationships  crucial  to  the  design  of  natural-star  and  synthetic-beacon  adaptive-optics 
systems  have  now  been  assembled.  A  comparative  evaluation  of  these  two  compensation  approaches  is 
most  easily  accomplished  through  a  calculation  of  the  achievable  long-exposure  Strehls.  Recall  that  the 
Strehl  for  uncompensated  turbulence  is 

[Strehl]^  «  — .  (4-21) 
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A  modest  improvement  to  this  number  can  be  realized  by  incorporating  a  tracking  system  to  remove  the 
tilt  component  of  turbulence,  thus  achieving  an  on-axis  irradiance  similar  to  the  short-exposure  Strehl. 
From  Yura  [14], 
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thus,  for  the  aperture  diameter  and  wavelengths  of  interest  this  relationship  predicts  a  Strehl  improvement 
of  approximately  50  percent.  To  obtain  more  substantia]  gains,  the  high  spatial  frequency  components 
must  also  be  corrected.  As  discussed  in  Section  2.4,  the  long-exposure  beam  profile  that  results  from 
partial  phase  compensation  is  predicted  to  have  a  residual  background  of  width  X/r{)  and  a  diffraction- 
limited  centra]  core  that  is  smeared  as  a  result  of  tilt  jitter.  If  the  ratio  between  the  uncorrected  back¬ 
ground  and  the  central  core  is  large,  the  Strehl  ratio  can  be  approximated  by  the  following  expression: 
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For  natural  beacons,  the  focus-anisoplanatism  term  in  the  numerator  will  be  zero;  in  addition,  tracking 
errors  in  the  denominator  can  be  ignored,  as  the  constraints  imposed  by  the  wavefront  sensor  require  a 
very  bright  reference  having  a  small  angular  displacement.  The  beacon-offset  term  in  the  numerator  will 
be  zero  for  synthetic-beacon  systems. 

The  strategy  at  this  point  will  be  to  develop  a  balanced  error  budget  based  on  Equation  (4-23).  In 
keeping  with  the  relatively  modest  goals  that  are  currently  being  discussed  within  the  astronomy  com¬ 
munity,  a  set  of  design  criteria  will  be  used  that  results  in  a  factor-of-10  ratio  between  the  central  core 
of  the  beam  profile  and  the  background  and  which  limits  the  width  of  the  central  core  to  twice  the 
diffraction  limit.  The  corrected  and  uncorrected  point-spread  functions  are  compared  in  Figure  24  for 
imaging  at  0.55  pm. 
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Figure  24.  Comparison  of  the  corrected  and  uncorrected  beam  profiles  for  a  design  criteria  that  restricts  the  tilt 
jitter  to  twice  the  diffraction-limited  beam  diameter  and  achieves  a  signal-to-background  ratio  of  10.  This  plot 
represents  the  point-spread  function  comparison  for  operation  at  0.55  pm. 


The  final  results  of  the  analysis  of  a  4-m  compensated-imaging  system  are  summarized  in  Table  3. 
A  total  of  six  configurations  are  included,  and  the  error  budget  for  each  can  be  found  in  the  uppermost 
section  of  this  chart.  The  corresponding  component  specifications  are  given  in  the  next  four  rows,  which 
show  the  required  number  of  actuators,  servo  bandwidth,  average  laser  power,  and  beacon  brightness.  It 
is  not  surprising  that  in  each  case  the  actuator  spacing  is  closely  matched  to  turbulence  coherence  length 
and  that  the  servo  bandwidth  (defined  as  1/lOr^)  is  approximately  equal  to  the  Greenwood  frequency. 
With  the  exception  of  the  requirement  on  the  laser  output,  all  of  the  numbers  listed  are  well  within  the 
current  state  of  the  art. 

The  bottom  section  of  Table  3  contains  a  quantitative  comparison  of  the  long-exposure  Strehl  and 
resolution,  as  defined  by  Equations  (2-49)  and  (2-50).  In  each  case  a  substantial  improvement  in  the 
effective  seeing  and  image  intensity  is  achieved.  The  motivation  for  pursuing  synthetic-beacon  constructs 
is  clearly  indicated  in  the  comparison  of  operational  utility  listed  in  the  last  row  of  this  summary.  As  a 
result  of  a  lower  irradiance  requirement  for  full-aperture  tracking  (corresponding  to  6  visual  magnitudes) 
and  a  larger  isoplanatic  angle  (by  a  factor  of  3.5),  the  probability  of  finding  a  reference  star  for  tilt 
stabilization  in  the  visible  is  over  3  orders  of  magnitude  greater  than  the  likelihood  of  finding  a  beacon 
source  for  phase  compensation.  Although  the  absolute  utility  at  0.55  |im  is  only  a  fraction  of  a  percent, 
because  this  number  scales  roughly  as  z qI^Sq,  a  slightly  more  benign  turbulence  model  would  result  in 
a  greatly  improved  estimate. 
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In  closing,  the  system  designer  need  not  be  constrained  to  precisely  tailor  the  adaptive-optics 
equipment  to  the  expected  turbulence  conditions  or  the  correction  wavelength.  Through  the  use  of  inter¬ 
changeable  reconstruction  operators  and  adjustable  timing  circuitry,  it  is  always  possible  to  modify  both 
the  effective  actuator  density  and  servo  bandwidth  to  achieve  maximum  effectiveness.  Thus,  the  short- 
wavelength  capability  of  the  system  can  steadily  evolve  as  component  technology  improves  without 
suffering  a  loss  of  performance  in  the  infrared. 


TABLE  3 

Summary  of  Compensation-System  Design  Specifications  and  Operational  Utility  for  a 
Factor-of-5  Improvement  in  the  Long-Exposure  Strehl  for  a  4-m  Telescope 


Natural-Star  System 

Laser  Guide-Star  System 

V  Band 
0.55  pm 

J  Band 
1.25  pm 

K  Band 
2.2  pm 

V  Band 
0.55  pm 

J  Band 
1.25  pm 

K  Band 
2.2  pm 

°ftgure  (rad2  of  Phase> 

4.2 

2.4 

1.3 

2.9 

1.2 

0.49 

(rad2  of  tilt) 

0 

0 

0 

1.2X10"14 

5.9  x10“14 

1.8X10"13 

Number  of  Actuators 

290 

85 

45 

250 

90 

65 

Servo  Bandwidth,  1/1 0xd  (Hz) 

23 

12 

9 

21 

12 

11 

Sodium-Beacon  Laser  Power  (W) 

__ 

— 

— 

47 

9 

6 

Beacon/Fiducial-Star  Magnitude 

8 

10 

11 

13 

16 

18 

Uncorrected  Strehl 

0.0016 

0.012 

0.044 

0.0016 

0.012 

0.044 

Corrected  Strehl 

0.016 

0.11 

0.32 

0.015 

0.082 

0.17 

Uncorrected  Resolution  (x  1.22  A/D) 

25 

9.2 

4.8 

25 

9.2 

4.8 

Corrected  Resolution  (x  1.22 A/D) 

2.6 

1.3 

1.0 

3.1 

2.0 

2.0 

Fractional  Sky  Coverage 

3X10-6 

7x10-5 

3X10-4 

0.003 

0.2 

1 
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5.  CONCLUSIONS 


The  analytical  relationships  developed  in  this  optimization  study  are  generally  consistent  with 
empirical  design  rules  that  have  evolved  over  nearly  three  decades  of  research  in  the  field  of  high-energy 
laser-beam  control.  For  imaging  applications,  the  requirements  for  phase  correction  can  be  relaxed  some¬ 
what,  but  to  achieve  a  significant  improvement  in  the  short-exposure  beam  profile,  the  actuator  spacing 
and  correction  bandwidth  of  an  adaptive-optics  system  must  be  closely  matched  to  the  corresponding 
turbulence  parameters. 

At  this  point  there  is  little  doubt  as  to  the  applicability  of  adaptive  optics  to  the  field  of  astronomy. 
Experiments  conducted  by  researchers  at  Lincoln  Laboratory  and  other  facilities  have  conclusively  dem¬ 
onstrated  the  efficacy  of  both  the  wavefront-compensation  technology  and  the  laser  guide-star  concept. 
Much  still  needs  to  be  done,  however,  before  adaptive-optics  technology  is  widely  accepted  as  a  tool  for 
astronomical  research.  Effectively  adapting  this  technology  to  astronomical  telescopes  will  require  an 
emphasis  on  reducing  overall  system  complexity  and  cost  while  maintaining  component  reliability.  The 
first  step  in  this  process  is  to  establish  acceptable  criteria  for  performance,  upon  which  a  suitable  balance 
between  capability  and  design  sophistication  can  be  achieved.  This  issue  has  been  the  principal  motiva¬ 
tion  in  the  analysis  just  presented. 
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APPENDIX  A— AN  ANALYTICAL  TREATMENT  OF  TURBULENCE 

EFFECTS 


Quantitative  estimates  of  turbulence  effects  have  traditionally  been  difficult  to  obtain  because  their 
analytic  treatment  involves  complex  integrals  and  fractional  powers  of  the  integration  variable.  As  a 
result,  most  investigators  have  resorted  to  numerical  integration  techniques  or  Monte  Carlo  methods  that 
use  a  large  ensemble  of  individual  turbulence  realizations.  Although  these  approaches  give  accurate 
results,  parametric  investigations  are  difficult  and  the  underlying  physics  is  often  lost. 

The  genera]  philosophy  of  the  Lincoln  approach  is  that  far-field  characteristics  (such  as  Strehl  and 
resolution)  can  be  derived  from  residual  phase  error  as  measured  in  the  pupil  plane  of  the  telescope.  If 
an  adaptive-optics  system  is  used,  residual  phase  error  refers  to  the  difference  between  the  actual  phase 
and  the  phase  correction  applied  by  the  deformable  mirror.  These  differences  can  arise  from  insufficient 
temporal  or  spatial  resolution  in  the  compensation  system  or  differences  between  the  beampath  of  the 
turbulence  beacon  and  the  object  under  investigation.  All  these  effects  can  be  studied  using  the  approach 
described  in  this  report. 

A.l  STATISTICAL  DESCRIPTION  OF  THE  PUPIL-PLANE  ELECTRIC  FIELD 

For  imaging  applications,  the  ultimate  interest  lies  in  the  point-spread  function  of  the  collection 
system  in  the  focal  plane,  but  for  turbulence  calculations  this  quantity  is  most  easily  inferred  from  the 
residual  phase  error  in  the  plane  of  the  receiver  aperture.  Although  amplitude  fluctuations  play  a  role  in 
the  distortion  process,  to  a  good  approximation  the  far-field  error  can  be  accounted  for  by  the  average 
phase  variance,  CT^,  of  the  electric  field  in  the  pupil  plane  of  the  telescope. 

A.1.1  The  Rytov  Approximation 

The  general  form  of  the  wave  equation  for  propagation  through  turbulence  cannot  be  solved  in 
closed  form,  but  through  the  judicious  application  of  approximations  an  approximate  form  can  be  devel¬ 
oped  that  is  amenable  to  solution.  The  usual  starting  point  is  a  form  introduced  by  Tatarski  (37]  that 
expresses  the  electric  field  as  a  function  of  the  three-dimensional  position  vector,  f ,  and  time,  /,  in  the 
presence  of  turbulence: 

E(r,l)  =  E1/2  exp[/^  •  r  -  i(Ot\  exp[*(F,/)  +  i0(r,t)]  (A-i ) 

Eq  Ej 

where  £0  =  unperturbed  radiation, 

£,  =s  perturbation  component, 

E^2  s  average  irradiance, 

(O  =  temporal  frequency  of  the  electromagnetic  field. 

Icq  =  wave  vector  of  magnitude  2 it/. A,  and 
A  =  wavelength  of  light. 


53 


The  real  parameters  X(rJ)  and  0(t\t)  define  the  log-amplitude  and  phase  components  of  a  multiplica¬ 
tive  perturbation  that  is  introduced  when  light  is  diffracted  as  a  result  of  changes  in  the  index  of  refraction 
along  the  propagation  path.  This  form,  known  as  the  Rytov  approximation,  is  valid  as  long  as  the 
perturbation  is  small.  From  Equation  (A-l)  it  follows  that 

£  +  =  Inf  1  +  — )  =  —  .  (A-2) 

{  Eo)  E0 

It  will  be  assumed  that  x  2nd  <f>  are  statistically  independent  random  variables  having  Gaussian  probability 
distributions  and  zero  mean. 

A  detailed  derivation  of  the  perturbation  equation  can  be  found  in  several  excellent  texts  [38]; 
therefore,  only  a  brief  outline  of  the  standard  treatment  of  this  problem  will  be  given.  As  indicated  in 
Figure  A-l,  light  passing  through  a  turbulent  medium  is  assumed  to  diffract  from  a  set  of  phase  sheets 
aligned  in  planes  perpendicular  to  the  propagation  axis.  If  the  z  axis  is  taken  as  the  propagation  direction 
and  the  effect  in  the  plane  of  the  origin  of  the  coordinate  system  at  time  t  =  0  is  measured,  the  perturbation 
component  of  the  electric  field  can  be  obtained  directly  from  the  equation  for  Fresnel  diffraction. 


Figure  A-l.  Propagation  of  light  through  a  single  phase  screen  at  range  z  '  to  an  aperture  of  diameter  D  -entered 
at  the  origin  of  the  coordinate  system. 
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The  function  rii(x\y\z')  represents  the  perturbation  component  of  the  refractive  index,  and  the  volume 
integration  is  performed  over  the  entire  space  between  the  origin  and  the  source.  By  combining  this 
expression  with  Equations  (A-l)  and  (A-2),  the  following  pair  of  equations  is  obtained: 


X(x,y,0)  =  -^°jdz']  dx'  J  dy’^ cos- 


^[(x-x'f  +  ly-yf 


2z' 


"I  (•*'./.  2'), 


(A-4) 


and 


2jti  J  J  z  2z 


nx(x',y',z') 


(A-5) 


The  integrals  over  x  and  y  are  recognized  as  the  convolution  of  a  filter  function  with  the  random 
variable  n,.  For  a  linear  transformation  of  this  form,  the  power  spectral  density  of  the  output  variable  is 
equal  to  the  product  of  the  squared  modulus  of  the  transfer  function  of  the  filter  and  the  spectral  density 
function  of  the  input  variable.  The  transfer  functions  of  the  filter  expressions  are  obtained  from  their 
Fourier  transforms 


(A-6) 


(A-7) 
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where  Kx  &2n/x  and  Xy  =  2Jljy  are  the  spatial  wavenumbers.  From  this  result  the  power  spectral 
density  functions  of  %  $  are  seen  to  be 


and 


0x^Kx,Ky)  =  k^  jdz'  sin1 


x2x  +  xy 
21* 


K„Ky-.z') 


(A-8) 


COS 


<1 


2k$ 


(A-9) 


Notice  that  the  spectral  density  function,  <Pnf)c‘I,K'y;z,)>  that  appears  in  the  last  two  equations  is  only 
defined  within  the  plane  in  ic  space.  The  z  parameter  is  explicitly  included  in  the  variable  list 
because  the  strength  of  the  refractive  index  fluctuations  are  typically  dependent  on  position  along  the 
propagation  path. 


A.1.2  Three-Dimensional  Spectral  Density  Function  Normalization 

At  this  point  in  the  derivation  it  becomes  necessary  to  develop  a  relationship  between  the  two- 
dimensional  spectral  density  function  and  its  three-dimensional  form,  0n^Kx,Ky,  in  which  the 

refractive-index  spectral  density  is  typically  presented.  In  this  representation  it  is  assumed  that  turbulence 
is  constant  locally,  and  its  average  properties  change  slowly  in  the  propagation  direction.  This  is  accom¬ 
plished  by  taking  a  careful  look  at  the  manner  in  which  the  z  integral  is  performed  and  the  effect  of 
correlation  along  the  z  axis. 

When  z-axis  correlation  effects  are  included  in  Equations  (A-4)  and  (A-5),  the  power  spectral 
density  of  an  output  variable,  a,  is  properly  written  as  the  convolution  of  the  filter  function, 
Fa^Kx,Ky;z'^F^Kx,  Ky;z with  the  cross-spectral  density  function,  Ky,z\z"},  of  the  refrac¬ 

tive  index 

oo  oo 

—  § dz  j dz  Fa^Kx,Ky,z}Fa{xx,Ky,Z  )  .  (A-10) 

0  0 

By  making  the  variable  substitutions  z  =  (z'  +  z")/2  and  Az-z’-z"-  this  expression  can  be  rewritten 


oo  Z 

0a{xx,  Ky)-jdzj dAzFa(Kx , Ky;z  +  Az/2) F*(kx, Ky;z-  Az/2) 0„( Kx, Ky;Az,z) 
o  -z 

oo  2  00 

*  jdz\Fa(Kx,Ky;z)\  jdAz0n(Kx,Ky;Az,z ) 
o  -oo  ’  (A- 11) 
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where  the  interchange  of  the  order  of  integration  and  the  modification  of  integration  limits  is  possible 
because  the  correlation  length  of  &n  is  short  with  respect  to  functional  changes  with  respect  to  e.  This 
equation  is  similar  in  form  to  Equations  (A-8)  and  (A-9),  and  the  cross-spectral  density  integral, 

oo 

j  dAz0n(Kx,Ky',Az,z)'  *s  seen  t0  b®  equivalent  to  the  two-dimensional  spectral  density  function.  The 
desired  result  is  obtained  by  noting  that  <P„  and  &„  are  Fourier-transform  pairs  with  respect  to  Az.  so 


1  °° 

!F^0n(Kx,Ky;Az,z)]  =  —  $dAz<i>n(Kx,Ky-,Az,z)cos(AzKz)  =  0„(Kx,Ky,Kz;z) 


(A-12) 


0n(Kx,Ky;z)=  jdAz<Pn(Kx,Ky;Az,z)  =  27t4>„(Kx,Ky,0-,z)  . 


(A-13) 


From  Equation  (A-13)  it  can  be  seen  that  the  expressions  for  log-amplitude  and  phase  can  be  rewritten 
in  a  form  that  is  compatible  with  a  three-dimensional  refractive-index  spectrum  by  inserting  the  multi¬ 
plicative  constant  2rt 


<Px(kx,  JCy)  =  2; 1$  J dz  sin2[  K-~y-z  <P„(kx,  vy,0;z) 
0  V  ZKq  ) 


(A- 14) 


0,p(kx,  Ky)  =  2nk£  j  dz  cos 


2f  Kx  +  Ky 


z  0n(Kx,Ky,O ;z) 


(A- 15) 


A.  1.3  The  Kolmogorov  Turbulence  Spectrum 

The  power  spectral  density  function  of  the  index  of  refraction  can  be  derived  from  well-established 
physical  laws  relating  to  the  turbulent  flow  of  air.  For  spatial  frequencies  within  a  region  known  as  the 
inertial  subrange,  Kolmogorov  [7]  showed  that  the  structure  function  is  proportional  to  the  2/3  power  of 
the  scalar  distance.  From  Kolmogorov’s  observation,  it  follows  that 


<Pn(ic,z)  =  0.033C%(z) 


(A- 16) 


2  2  2  2  2 

where  V  =  Kx  +  Ky  +  Kz ,  and  C„(z)  is  the  strength  of  the  turbulence  fluctuations  along  the  propaga¬ 
tion  path.  As  shown  in  Figure  A-2,  at  large  spatial  frequencies  the  spectrum  decays  more  rapidly  than 
predicted  by  Equation  (A- 16)  due  to  the  dissipation  of  small  eddies  by  viscous  forces.  Furthermore,  at 
small  spatial  frequencies  the  function  does  not  become  infinite  but  approaches  a  constant  value.  A  more 
precise  representation  is  given  by  the  von  K  arm  In  spectrum 
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Figure  A-2.  The  von  Karmdn  turbulence  spectrum.  The  Kolmogorov  spectrum  is  valid  within  the  inertial  subrange. 


4>n{X'Z) 


0.033C2(z) 


/ 

exp 

V 


(A-17) 


where  L0  =  2k/ko  represents  the  outer  scale  of  the  inertial  subrange  and  K",  is  related  to  the  inner  scale. 
This  particular  form  for  the  inner  scale  effect  was  suggested  by  Tatarski;  other  models,  such  as  the  Hill 
spectrum  [39,40],  can  also  be  used.  For  most  ground-to-space  applications,  however,  sufficient  accuracy 
is  achieved  with  the  simpler  Kolmogorov  spectrum. 


A.1.4  Log-Amplitude  Variance  and  Phase  Variance  of  Uncorrected  Turbulence 

By  combining  Equations  (A*  14),  (A- 15),  and  (A- 16),  expressions  for  the  net  variance  of  the  pupil- 
plane  phase  and  log  amplitude  fluctuations  due  to  turbulence  are  obtained; 


oo  oo 


0.207 /to2  jdzj  dKx  J  dKyC2(z)K~u/3 sin: 


0  “ o 


2fco 


(A- 18) 


crj  =  0. 207  ^JdzJ  dKx  J  dKy  C2 (z)  K~l  1/3  cos2 


^  K2z^ 


0  — ° 


2ko 


(A- 19) 
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where  k2  =  K2  +  ■  (In  all  subsequent  calculations  the  Jf-space  integration  will  be  confined  to  the  k*.  kv 

plane.)  These  two  equations  represent  the  standard  representation  of  the  total  variance  due  to  uncorrected 
turbulence.  The  following  section  will  show  how  this  formalism  can  be  extended  to  address  a  much 
broader  range  of  turbulence  problems. 


A.2  A  GENERALIZED  APPROACH  TO  TURBULENCE  CALCULATIONS 


The  essentia]  contribution  of  this  work  derives  from  two  important  observations.  The  first  discovery 
is  that  a  wide  variety  of  turbulence-related  problems  can  be  treated  analytically  by  inserting  an  appro¬ 
priate  multiplicative  filter  function,  F(k\z),  into  Equations  (A- 18)  and  (A- 19).  When  a  single  beam  is 
involved,  the  general  expression  for  phase  variance  has  the  form 


<rj  =0.207*o2  ldzC2n(z)  \  dKK~m  F{k,z)  ,  (A.20) 

0 

where  the  approximation  cos2(jc)  =  1  can  be  made  for  most  problems  in  which  the  wave  propagates  from 
space  to  ground  without  introducing  significant  error. 

The  integrand  of  the  triple  integral  indicated  in  the  previous  expression  will  usually  have  a  very 
complicated  structure  that  includes  fractional  powers  of  k  and  Bessel  functions  of  k  and  z.  Our  second 
observation  is  that  the  k  integration  can  be  treated  as  a  Mellin  transform,  tables  of  which  can  be  found 
in  standard  mathematical  handbooks  [41].  The  Mellin  transform  of  g(x)  is  defined  as 


G(s)  =  M[g(x)]  =  jdxxs~'g(x)  .  (A-2I) 

o 

Thus  an  integral  over  k  that  contains  the  term  K?  is  equal  to  the  Mellin  transform  of  the  remaining 
integrand  evaluated  for  s  =  p+ 1 ; 

oo 

jdKKpg(K)  =  G(s)\s=p+l  .  (A-22) 

0 

The  application  of  these  computational  techniques  is  best  illustrated  by  a  few  simple  examples.  For 
simplicity,  only  the  phase  component  of  the  electric  field  will  be  explicitly  treated  in  the  next  section. 
However,  the  extension  of  these  methods  to  calculations  involving  the  log-amplitude  component  is 
straightforward. 

A.2.1  Zernike  Components  of  the  Phase  Variance 

The  Zemike  modes  of  a  circular  aperture  are  incorporated  in  virtually  all  calculations  of  turbulence 
distortion.  The  expression  for  the  variance  of  the  i ^  Zemike  component  of  the  phase  function  over  an 
aperture  of  diameter  D  is 
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In  D/2 

j  d6  JrfppZJ(p10)0(p,O) 


o  o 


(A-23) 


where  Z,(p,0)  is  the  Ith  Zemike  polynomial.  Because  the  statistical  properties  of  the  phase  are  defined 
by  its  power  spectral  density  function,  this  integral  is  more  easily  evaluated  in  the  Fourier  domain.  Fourier 
transforms  of  the  Zemike  polynomials  evaluated  for  an  aperture  radius  of  unity  are  found  in  Noll  (42): 


•  jL  Ji.  i 

Qi{ky<p)  =  ~  |  d9jdppZi{p,$)cxp^—2nik-pj  .  (A-24) 

*o  o 

These  functions  can  be  adapted  to  the  notation  used  in  this  report  by  making  the  transformation 
k  -  kD/4k,  so  that 


2  n 


{(7^  =  0.207*0  jdzCl(z)  J  d(p\  dKK  8/3|Q,(!cD/4;r,<p)| 


(A-25) 


As  these  components  are  orthogonal,  the  phase  variances  associated  with  different  Zemike  modes  will 
add  incoherently. 

Piston-Removed  Phase  Error  for  Uncorrected  Turbulence.  The  Kolmogorov  spectrum  is  not  only 
infinite  at  the  origin,  its  total  phase  variance  is  also  infinite.  This  result  is  non-physical  and  can  be 
remedied  by  applying  a  finite  outer  scale;  however,  the  more  appropriate  approach  is  to  recognize  that 
an  average  phase  shift  over  the  collection  aperture  will  produce  no  net  distortion  in  the  far  field  and 
should  therefore  be  ignored.  The  aperture-average  phase  is  referred  to  as  “piston”  by  the  adaptive-optics 
community,  and  its  value  is  derived  from  the  product  of  the  first  Zemike  polynomial  and  the  pupil-plane 
phase.  The  filter  function  for  piston  removal  over  an  aperture  of  diameter  D  is  found  to  be 


F(k)  =  1  -  j  (2i(k  D/4/t)|2  =  1  - 


2Jx{kd:i) 
kD/2 


(A-26) 


The  K  integration  specified  by  Equation  (A-20)  is  performed  in  the  following  manner: 


c.*--n/3 


dKK 


2J\{kDI2) 
kD/2 


l2l 


=  2nj  dKK 


-8/3  I 


2Jx(kD’2) 
k  D/2 


121 


D2 

J 

(A-27) 


Lr=-1 1/3 
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The  last  expression  is  similar  in  form  to  a  Mellin  transform  given  in  Appendix  B: 


i  r[i  +  ^,2]r[i,2-5-2] 

2\n  rp-s:2]r\l-s  2] 

By  applying  the  transform  identities 


M[-g(x)]  m -M[g{xj\  and  <M[g(ax)]  s  a  s  M[g(x)\ 
the  following  relationship  is  obtained: 


(A-l>8> 


(A -29) 


32  n 
D2 


(kD/2Y 


-J?(kD/2) 


32 1  r[- 5  6]r[l4-6] 


r=-l  1/3 


4  J 

=  2.11  D5'3 

Combining  this  expression  with  Equation  (A-20)  yields 


D2  \  2  )  2\  r  r[32  6]r[l7  6] 


(A-30) 


=  0.437 £>5'3  f dz  C2 (z)  ,  (A-31) 

o 

The  final  integration  over  z  cannot  be  completed  without  specifying  the  turbulence  profile;  how¬ 
ever,  this  calculation  is  easily  performed  for  the  standard  turbulence  models.  It  is  common  to  characterize 
a  turbulence  model  by  its  moments  fln  =  J  dhhn  Cn  (h),  which  are  computed  at  zenith.  For  an  arbitrary' 
zenith  angle,  z. 


jdzzn  C2(z)  s  secn+,(£)  jdhhn  C2(h)  =  secn+1(g)^„  .  (A-32) 

o  o 

With  the  application  of  thi<=  identity, 

<rj  =  0  437 D5  3  sec(C)^0  ,  (A-33) 

which  is  the  standard  result.  Although  this  particular  expression  can  be  easily  obtained  without  the  use 
of  Mellin  transforms,  the  benefit  of  this  approach  becomes  obvious  when  more  complicated  filter  func¬ 
tions  are  involved. 

Figure  and  Tilt  Calculations  for  Uncorrtcted  Turbulence.  Most  adaptive-optics  systems  incorpo¬ 
rate  separate  devices  for  tracking  and  high-spatial-frequency  compensation;  a  typical  cor  figuration  is 
shown  in  Figure  A-3.  For  this  reason  it  is  convenient  to  develop  individual  expressions  for  each  of  these 


61 


turbulence  components.  The  filter  function  for  tilt  is  derived  from  the  Fourier  transforms  of  the  second 
and  third  Zemike  polynomials; 


(  kD;2  J  (,  K  0/2 


\2' 

'AJ2(kD'2) 

/ 

kD  2 

(A-34) 


and  the  associated  Mellin  transform  is 


sv/jyfW] 


1  r[2  +  s/2]r[l/2-s/2l 

2v'7r  r[3-s'2]r[i-s:2] 


(A-35) 


The  tilt  component  of  turbulence  is  usual  y  given  in  units  of  rad2  of  single-axis  tilt  angle,  and  the 
conversion  from  rms  phase  to  rms  angle  involves  the  multiplicative  constant  8 /(k^j2.  Combining  this 
constant  with  Equation  (A-35)  and  the  variable  change  indicated  in  Equation  (A-29)  results  in 


crjju  =  3.04  D  13  sec(£)// 0  (rad2  of  single-axis  tilt) . 


(A-36) 


The  piston-  and  tilt-removed  phase,  also  known  as  “figure,”  represents  the  high-spatial-frequency 
turbulence  that  must  be  corrected  by  the  deformable  mirror.  The  filter  function  for  figure  variance  is  the 
difference  between  Equations  (A-26)  and  (A-34); 


^figure  (*") 


’U^kD/I) 

2 

~  A  J2(k  D/2)12 

kD/2 

kD/2 

By  applying  the  identity 

g}  (x)  +  g2(x)]  s  94[gx (x)]  +  94[g2{x)\  , 


(A-37) 


(A-38) 


the  following  result  is  obtained; 


94 


•/iW- 


+  4  J\(x) 


1  r[l  +  j/2]r[l/2-j/2]  2  r[2  +  s:2]r[\'2-s/2 ] 

2 r[2-s/2]r[l-V2]  *  \rn  r[3~ s/2]  f[l  -s  2]  -  (A-39) 
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Figure  A-3.  Schematic  presentation  of  a  closed-loop  adaptive-optics  system.  The  first  derivative  of  the  phase 
difference  between  the  incoming  wavefront  and  the  surface  figure  imposed  by  the  tilt  and  deformable  mirrors  is 
measured  by  the  phase  sensor.  The  residual  phase  error  is  computed  by  the  reconstructor  and  subsequently  applied 
to  the  two  active  optical  elements. 


The  final  expression  for  figure  variance  is  derived  by  following  the  steps  outlined  in  the  previous 
examples: 


<^e»«  -  0.207 *o2  }*c,’k)  0?*)  M 


=  0.0567^  D5/3sec(£)jt0 


■/|2(kD/2)- 


[kD!  2) 


+  4jf(fcD/2) 


li=-l  1/3 


(A-40) 
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A.2.2  Phase  Variance  due  to  Anisopianatic  Effects 

Most  problems  of  interest  for  compensated-imaging  applications  involve  a  pair  of  beams,  one 
which  originates  from  the  beacon  that  is  viewed  by  the  wavefront  sensor  and  a  second  that  passes  through 
the  correction  system  and  is  imaged  onto  the  focal  plane.  Anisopianatic  effects  occur  whenever  these  two 
beams  follow  different  paths  through  the  atmosphere.  Figure  A-4  shows  a  simple  geometry  in  which  the 
rays  impinging  on  a  point  in  the  pupil  plane  of  the  receiver  pass  through  points  in  a  single  phase  screen 
that  differ  by  the  vector  displacement  d .  If  an  adaptive-optics  system  is  used  to  apply  a  correction  based 
on  the  beacon  information,  the  magnitude  of  the  resulting  error  at  point  p  will  be  proportional  to 

|n,(p;2)-nj(p  +  d(p;z);z)|2  ^ 

where  it  is  noted  that  d  will  be  a  function  of  both  aperture  position  and  the  range  to  the  phase  screen. 
In  K  space  this  displacement  corresponds  to  a  phase  shift  of  K-d,  so  that  the  corresponding  error 
expression  would  be 

1 1  —  exp[r( K  ■  d{jp\  z))]  |  |7[n|(p;2)]|2=2[l-cos(Y-rf(p;z))]^n(v;-)  _  (A-41) 

The  term  2|l  -  cos(fC- d(p;z))j  can  be  viewed  as  a  filter  function  associated  with  the  displacement  of 
the  two  rays.  To  obtain  a  function  representative  of  the  aperture-average  error,  the  following  integral  must 
be  performed: 

4  ln  m 

F(^)  =  -p- J  de  J dp  p  [2  -  2  cos(  k  ■  d(p,  0;  z))]  .  (A-42) 

0  0 

For  some  viewing  geometries  it  may  be  possible  to  develop  the  filter  function  ass  "iated  with  a 
single  component  of  an  anisopianatic  effect  by  taking  the  product  of  Equation  (A-42)  and  ti.o  appropriate 
Zemike-function  filter.  For  the  more  general  case  illustrated  in  Figure  A-5,  however,  the  problem  in¬ 
volves  two  beams  that  may  have  different  diameters,  intensity  distributions,  spatial  offsets,  and  may 
originate  from  sources  placed  at  different  altitudes.  The  pupil-plane  error  variance  results  from  the 
difference  function 


^  =  0.207/^  jdzj  dp  j  dK&n(ic,z) 

,  [cos[^(ri,  *;  z)]1  ,  _Jcos[P2(y2,  v;z)]l 

X  C>M  •  for  M  ~A(^z)G2(y2k)  .  r_  (  , 

|_sin[*i(yi,  K ;  z)j  J  [stn[p2(y2.  -)]  J 

ic°sh(ri,M]]  Jcos[/>2(y2,*;z)]1 

[sinjp^y,,  jc;  z)J  J  [sin[A(y2,  z)]  J 


,  (A-43) 
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J073t  *3» 

OBJECT  BEACON 


Figure  A-4.  Anisoplanatic  error  occurs  when  the  beacon  rays  that  impinge  on  a  point  in  the  collection  aperture 
follow  a  different  path  than  the  rays  from  the  object  under  observation.  Within  a  phase  screen  at  range  z the  rvi  o 
rays  pass  through  points  that  are  separated  by  the  vector  displacement  d. 


which  defines  both  the  phase  and  log-amplitude  error.  For  space-to-ground  propagation,  the  following 
parametric  definitions  apply: 

0„(k;z)  =  turbulence  spectrum, 

A(k,p,z)  -  amplitude  ratio  of  the  second  beam  relative  to  the  first, 

P(y,  K-,  z)  =  yK2 z/2kQ  =  propagation  parameter, 

7  =  (L-z)/L  for  a  source  positioned  at  range  L,  and 

G(yic)  =  jdr' g(r')e\p[iyic-r']  =  aperture  filter  function. 

It  is  important  to  recognize  that  the  character  of  the  the  amplitude  ratio  will  depend  on  the  specification 
of  the  aperture  filter  function.  When  G1  and  G2  are  unity,  A(ic,p,z)  describes  the  ratio  between  each  pair 
of  phase  points  at  range  z  that  lie  along  rays  intersecting  the  collection  aperture  at  point  p  ;  in  this  case 
an  integration  over  p  is  required  to  obtain  the  average  variance.  If  filter  functions  representing  an 
aperture-average  Zemike  component  (such  as  piston  or  tilt)  are  used,  however,  the  amplitude  ratio  is  then 
computed  only  for  the  beam  centroids,  in  which  case  v4(  will  be  a  function  of  z  but  independent  of  p. 
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With  the  introduction  of  Equation  (A-43),  two  quantities  can  be  computed  that  are  of  great  interest 
to  the  astronomical  community.  The  first  relates  to  the  issue  of  tracking  stability,  and  the  second  defines 
the  residual  figure  error  that  can  be  expected  with  the  application  of  laser  guide  stars. 
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BEAM  #1  BEAM  #2 


Figure  A-5.  Many  turbulence  calculations  involve  two  separate  beams.  For  adaptive -optics  correction,  one  beam 
may  represent  the  light  from  a  beacon,  whereas  the  second  beam  is  the  light  received  from  a  celestial  object. 


TiU  Component  of  Offset  Anisoplanatism.  For  situations  in  which  the  object  under  observation  is 
too  dim  to  track  directly,  a  bright  neighboring  star  might  be  used  as  the  reference  beacon.  However,  this 
geometry  will  give  rise  to  a  differential  tilt  due  to  the  angular  displacement,  #,  between  the  two  objects. 

As  indicated  in  previous  discussions,  the  aperture  filter  function  for  single-axis  tilt  is 


^tilt  ( *)  ~ 


(  \ 

8 

4J2(kD/2) 

l(*oO)2J 

kD/2 

(A-44) 


The  effect  of  an  angular  offset  of  0  is  introduced  into  the  variance  expression  through  the  amplitude 
function 

A(£-;z)  =  exp[iv- t?z]  ,  (A-45) 

which  represents  a  phase  shift  due  to  a  beam-centroid  displacement  that  increases  linearly  with  range. 
Combining  the  last  two  relationships  with  equation  (A-43)  produces  the  following  expression  for  tilt 
variance  for  Kolmogorov  turbulence: 


66 


w  4*S\  - 

<4  =  0.207  jdzCt(z)  J  dq>\dKK 


2  it 

J‘ 

0  0 


,-8/3 


V 


4J2(kD/2) 
{k^Dfl  KD/2  J 


-|2 

2[l  ~cos{k'i?^cos(<p)}]. 


(A-46) 


The  integral  over  <p  is  accomplished  through  the  application  of  the  Bessel  function  identity 


Jn(x)  = 


which  yields 


(A-47) 


<4  =  1.30 kZ\dzC2n{z)  J dKK~m 


8 

4J2(kD/2) 

kD/2 

[2- 2-/0(ict9z)] .  (A-48) 


Although  the  K  integration  indicated  in  this  equation  can  be  performed  directly  using  Mellin  transforms 
that  result  in  an  expression  containing  generalized  hypergeometric  functions,  for  small  displacement 
angles  the  correct  result  can  be  obtained  by  first  expanding  [l  -  J0(vt?z)]  in  a  Taylor’s  series.  The 
integral  for  the  first  term  of  the  expansion  is  then  obtained  using  the  Mellin  transform  for  J2(x)  given 
in  Equation  (A-35)  for  s  =  -5/3 

<4  *  5.34s ec3(^)D_7/3//2  t?2  (rad2  of  single-axis  tilt) .  (A-49) 

This  approximation  is  reasonably  accurate  for  d  <  (D/4  X  104  m).  When  normalized  to  the  diffraction- 
limited  beam  diameter,  the  beam  jitter  is  seen  to  be  proportional  to  sec3,2(£)  A~’  D~1/6  fi22  &■ 

Figure  Component  of  Focal  Anisoplanatism.  Laser  guide  stars  are  required  to  perform  high- 
spatial-frequency  compensation  on  objects  that  are  too  dim  to  provide  phase  information  to  the  wavefront 
sensor  of  an  adaptive-optics  system.  This  compensation  geometry  results  in  an  error  known  as  focal 
anisoplanatism,  and  is  one  of  the  principal  errors  associated  with  guide-star  imaging.  As  indicated  in 
Figure  A-6,  the  turbulence  above  the  beacon  is  unsampled  and  therefore  completely  uncorrected.  The 
residual  figure  error  associated  with  this  high-altitude  component  is  taken  from  Equation  (A-40): 


=  0.0567*2  D!-Vc 

where  the  notation 


(A-50) 


^[H),\dhh-Cl{h) 

H 


(A-51) 


has  been  introduced  to  represent  the  upper  fractional  moment  computed  from  the  beacon  altitude.  H,  to 
the  top  of  the  atmosphere. 
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Figure  A-6.  The  vertical  rays  represent  light  originating  from  a  distant  source  that  accumulate  phase  error  in 
traveling  through  the  atmosphere  to  the  telescope.  Radiation  from  the  laser  guide-star  follows  a  slightly  different 
path  to  the  telescope  and  is  therefore  unable  to  sample  turbulence  distortions  above  the  beacon  altitude  or  to 
precisely  sample  the  turbulence  that  lies  below  the  beacon. 


Below  the  beacon,  the  object  and  beacon  rays  that  strike  a  point  p  in  the  collection  aperture  are 
displaced  from  one  another  by  a  distance 


rf(p;z)  =  ^  (A-52) 

at  a  range  z  from  the  pupil  plane.  Therefore,  from  Equation  (A-43).  the  total  phase  error  due  to  the 
turbulence  below  the  beacon  is 


L  lit  °o  .  In  0/2  r  . 

=0.207^ \dzCl(z)  J  d<pjdK  K~m J  d6  jdpp  2-2cos  -^cos(S) ) 

n  on  n  n  L  k  L  )_ 


0  0 


0  0 


(A-53) 
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The  integrals  over  6  and  p  can  be  performed  by  applying  the  Bessel  function  identities 

1  r  vn+l 

Jn(x)sT~  j</0cos[xcos(0)  +  u0]  and  J  dxxn+lJn(ax)  = - Jn+l (ax)  ,  (A-54) 

2lt  o  a 

'A  ith  the  result 


L  2n  oo 

<rj  =  0.207^  jdzCfrz)  J  d(p\dK  K 


-8/3 


2-4 


J{{kDz2L j 
kDz/2L 


0  0  0 

The  integration  over  k  is  accomplished  by  applying  the  Mellin  transform 


(A-55) 


'■  w-f 


yi  rfi/2+y2l 

r[3/2-i/2] 


(A- 56) 


for  s  =  -8/3.  In  the  final  expression,  all  the  integrals  have  been  performed  analytically  except  for  the 
fractional  moment  computed  between  the  receiver  and  the  beacon  altitude 


«}-0 ' 


(A-57) 


Equation  (A-57)  represents  the  total  error  due  to  the  turbulence  below  the  beacon.  To  compute  the 
figure  component  of  this  error,  the  piston  and  tilt  terms  must  be  removed.  If  the  beacon  is  placed  over 
the  center  of  the  receiving  aperture,  the  amplitude  ratio  A(ic;z)  will  be  unity,  but  the  values  of  y  for  the 
object  and  beacon  beams  will  be  1  and  l-z/L,  respectively.  Thus,  from  the  aperture  filter  functions  given 
earlier  in  Equations  (A-26)  and  (A-34), 


crj  =  1.30^J<feC„2(z)  jdK  K~m4v: 


Jv(kD/2)  Jv(kD(1-z/L)/2 ) 
kD/2  kD(1-z/L)/2 


(A-58) 


where  the  value  of  v  is  1  for  piston  and  2  for  tilt.  A  Taylor’s  series  expansion  of  the  second  Bessel 
function  and  its  denominator  about  the  point  kD/2  yields  an  expression  for  which  the  leading  term  is 
quadratic  in  the  expansion  parameter  z/L  ; 


crj  =  1.30 kt\dzC2n(z)]dK  ,c-8/3  4  v2  (j-) 


Jv(kD/2)  +  (kD;2)J'v(kD/2) 

o  0  ~  KD/2 

This  can  be  further  simplified  using  the  identity  nj„(x)~ xJ'n(x)  =  xJn+i(x)  to  give 

<rj  =1.30*o2  jdzCt(z)  ]<*  k~8'3  4  v2  f—l  [-/v+,(vD/2)  +  (v-l)-v^^ 
o  o  L 


(A-59) 


.  (A-60) 
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The  ic  integration  is  performed  using  the  Mellin  transform  identities  found  in  Appendix  B,  from  which 
the  following  two  quantities  are  derived: 


and 


0^  *  0.0833D5  3  kfi  sec(£)  for  the  piston  component, 

H 


0$  =  0.354  D5  3  kg  sec(£)^2-y^-  for  the  tilt  component. 

H 


(A-61) 


(A-62) 


Combining  these  results  with  Equation  (A-57)  yields  the  figure  variance  due  to  turbulence  that  is  improp¬ 
erly  sampled  below  a  laser  guide  star 


flower  =  £)5/3^Sec(C) 


0.500— -  0.437 


H 


5/3 


H z 


(A-63) 


The  net  figure  error  that  would  result  from  a  beacon  placed  at  an  altitude  H  over  the  center  of  the 
collection  aperture  is  the  sum  of  Equations  (A-50)  and  (A-63); 


^figure  ~  £)5'  3  ^0  sec(£) 


0. 0567 jUq  ( //)  +  0. 500 


Ps/siH) 


H 


5/3 


-0.437 


dim 


H1 


.  (A-64) 
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APPENDIX  B— RELATIONSHIPS  APPLIED  IN  THE  ANALYTICAL 
TREATMENT  OF  TURBULENCE 

B.l  PROPERTIES  OF  THE  MELLIN  TRANSFORM 

The  Mellin  transform  C(s)  of  the  function  g(x)  is  defined  as 

oo 

G(s)3iW[g(*)]3  jdxx^gix)  .  (B-l) 

o 

The  structure  of  this  relationship  yields  the  following  set  of  identities,  which  are  frequently  invoked  in 
the  application  of  Mellin  transforms  to  turbulence  problems: 

M[-g(x)]  s  -M[g(x)]  ,  (B-2) 

M[g{ax)]  =  a~s  !M[g(x)]  ,  (B-3) 

and 

•  (B-4) 

It  can  also  be  shown  that  the  convolution  of  two  functions  can  be  derived  from  the  product  of  their 
transforms; 


J— gi(x)g2(x/y) 

*  V 


=  ^[g\{x)\M[g2{x)\ 


(B-5) 


The  principal  benefit  of  the  Mellin  transform  for  turbulence  calculations  derives  from  the  ready 
availability  of  standard  transform  tables  that  encompass  a  wide  range  of  interesting  problems;  see  Table 
B-l.  The  solutions  typically  involve  the  ratio  of  the  product  of  gamma  functions,  and  as  a  result  the 
following  notation  is  adopted: 


r 


CXj,...,CCm 


rm-nPn)  • 


(B-6) 


Because  the  only  singularities  of  this  ratio  occur  when  the  arguments  of  the  gamma  functions  in  the 
numerator  are  negative,  pole-residue  integration  can  be  used  to  obtain  a  power-series  expansion  for  all 
expressions  of  this  form. 
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TABLE  B-1 


Partial  List  of  Meilin  Transforms 


g(x) 

cxp(-x) 

sin(jr) 

cos(jc) 

Jv(x) 


J2v(x) 


2vv! 


Jl(x)-  x 


2  v 


M[g(x)]s  jdxx$~lg(x) 


2  s~\hr 


2  s~ly[Jcr 


2S~  lr 


r[s] 

1/2  +  s/2 

L  1  ~  s/2  J 

s/2 

1/2  -  s/2 

s/2+v/2 

V/2+1-J/2 


1 


s/2+v,l/2-s/2' 
v  +  l-s/2,l-s/2 


Ju(x)Jk{x)  21  Jt 


Jv(x)- 


2  Ik 

s/2  +  i~,^-s/2,l-s/2 


\-sl2  +  t~A-s/2  +  ^l~s/2  +  ~£j 


2  s~'r 


s/2+v/2 

v/2+1-5/2 


2zv(v!) 

Ju(x)Jx(x) 

vjU+* 


2yfx 


r 


si  2+  v,  1/2-5/ 2 
[v +  1-5/ 2, 1-5/2 


2^+a^!A! 


2yfn 


r 


5/2  +  ^~4-5/2,l-J/2 

1  -5/ 2  +  ^,  1-5/2  +  ^,  1-5/2  +  ^- 


2  J 


5  >  0 

-1  <5<  1 

0<5<  1 

-v<  5  <  3  /  2 

-2  v<  5  <  1 

-J/+A<5<  1 

-V-2  <  5  <  -V 

-2v-4  <  5  <  -2v 

--/i-A-4  <  5  <  -//-A 


B.2  APERTURE  FILTER  FUNCTIONS 

The  filter  functions  used  to  obtain  the  Zemike  modes  for  a  circular  aperture  are  listed  below.  Also 
given  are  the  filters  for  the  second-moment  quantities  obtained  with  a  finite  receive  aperture  and  a  finite- 
size  source. 
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B 2.1  Zernike  Modes 


The  general  formula  for  the  filter  function  associated  with  a  Zemike  component  on  a  circular 
aperture  of  diameter  D  derives  from  Noll  [42]; 


Fn.A*) 

FoM)  j 

The  piston  phase-variance  filter  function  is 
T  o/tVn/V'’2 

F(ic)  = 


2Jn+](KD/2) 
kD  12 


^2 


2cos2(m<p) 
2sin2(»ip) 
1  (m  =  0) 


(B-7) 


2Jj(kD/2) 
kD  12 


(B-8) 


which  must  be  multiplied  by  (l/k^)2  to  obtain  the  mean-squared  value.  The  function  used  to  determine 
the  x,  y,  and  total  phase  variance  associated  with  aperture  tilt  is 


Fv(k) 

F(*)l 


4J2(kD /  2) 
kDI  2 


cos 

•_  2 


'{(p) 


sur(tp) 


(B-9) 


from  which  the  variance  for  the  single  axis  tilt  angle  can  be  obtained  by  multiplying  by  . 


In  some  problems  the  variance  must  be  calculated  with  some  of  the  Zemike  modes  removed.  This 
is  easily  accomplished  by  subtracting  the  phase  variance  due  to  these  components  from  the  total  variance. 
For  instance,  the  filter  function  to  remove  piston  and  tilt  variance  from  a  single  wave  is 


p(f)=i- 


2Jx{kD  !  2) 

2 

4J2(kD /  2)~ 

kD!2 

KD(  2 

(B-JO) 


B.2.2  Gradient  Tilt 

The  gradient,  or  G-tilt,  is  equal  to  the  average  phase  gradient  over  the  receiver  aperture.  The  filter 
function  for  G-tilt  variance  over  a  circular  aperture  is 
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M*) 

*t*). 


J](kD/  2) 


cosz((p) 

sin2(<p) 

1 


(B-U) 


The  single-axis  tilt-angle  variance  is  obtained  from  the  above  by  multiplying  by  Mk^D), 


B.23  Filter  Functions  For  Finite-Size  Apertures  And  Distributed  Sources 


The  coherent  average  of  a  quantity  from  a  point  source  at  a  range  L  over  a  receive  aperture  is  equal 
to  the  piston  of  the  quantity,  which  is  given  by 


F(k)  = 


Jx(kD{\~z!  L)n) 


kD(\-z/L)I2 

where  z  is  the  range  variable.  If  the  point  source  is  at  a  very  long  range,  this  reduces  to 


F(*)  = 


Jx(kD/2) 

kDI2 


The  filter  function  for  an  incoherent  source  of  diameter  Ds  is 


F(k)  = 


Jx{kDsz!2L) 

kD'Z/2L 


(B-12) 


(B-13) 


(B- 14) 


The  filter  function  for  an  incoherent  source  of  diameter,  Ds,  received  coherently  by  a  receive  aperture  of 
diameter  D,  is  the  product  of  Equations  (B-12)  and  (B-14); 


F(ic)  = 


\jx{kD{\-zIL)I2): 

1\j,{kDsz!2L) 

kD{\-z/L)/2 

kDsz/2L 

(B- 15) 


B.2.4  Filter  Functions  For  Focal  Anisoplanatism 


The  filter  function  for  the  difference  of  two  waves  can  be  found  by  integrating  the  phase  over  the 
source  and  the  aperture.  The  function  for  focal  anisoplanatism  with  a  uniform  circular  source  offset  b 
from  boresight  is 


rm-1  Z2J'  (KDrJ2L)2Jl(KDzl2L)J  f  Kbz\ 
1  '  kD'ZI  2L  kDz/2L  °\  L  ) 


(Kbz\ 

\Jx{kDsz!2L)' 

l  L  ) 

kDsz/2L 

-a 


(B-16) 
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The  filter  function  for  focal  anisoplanatism  with  a  uniform  circular  source  on  boresight  is 


F(k)  =  \ 


22J1(kDsz/2L)2J1(kDz/2L)  < 


kDsz/2L 


kDz/2L 


2  J1(kDsz/2L) 
kD.z  /  2L 


(B-17) 


If  the  source  diameter  is  zero,  then  one  obtains  the  filter  function  for  focal  anisoplanatism  with  circular 
symmetry  for  a  point  source  with  no  offset 


f(?)  =  2 


|  ^J,(kDzI2L) 
kDzI2L 


(B- 18) 


The  filter  function  used  to  find  the  piston  or  tilt  component  of  focal  anisoplanatism  with  a  uniform 
circular  source  offset  b  from  boresight  is 


*(*) 


2vJv(kD/2) 

kD/2 


-exp 


ikb 


2vJv{kD[\-z!  L}!2)  2J](kDsz/2L ) 
(kD[\-z/ L]/2)  ( kDsz/2L ) 


(B- 19) 


B 3  ERROR  PROPAGATOR  FUNCTION  FOR  SECTION  STITCHING 

The  phase  reconstruction  process  can  be  described  by  the  matrix  product 

0  =  |(A,A)“lAjg  =  B|  ,  (B-20) 

in  which  g  is  a  vector  composed  the  j:  and  y  gradient  measurements,  0  is  the  phase  vector,  and  A  is  the 
matrix  is  defined  as  follows: 


g  =  A±  .  (B-21 ) 

If  the  gradient  errors  are  uncorrelated,  the  error  propagator  (defined  as  the  relationship  between  the 
gradient  error  variance  and  the  error  variance  of  the  reconstructed  phase)  is  given  by 


Error  Propagator  = 


N, 


-Trace 


(B-22) 


where  N0  is  the  number  of  elements  in  the  phase  array.  The  error  propagator  is  approximately  unity  for 
a  two-dimensional  system  composed  of  100  phase  points  [16]. 

The  propagation  of  low  spatial-frequency  errors  due  to  section-tilt  fluctuations  in  a  multiple-beacon 
system  can  be  viewed  as  a  reconstruction  process  in  which  the  section  tilts  are  viewed  as  local  gradients 
having  a  non-zero  correlation.  In  this  case  the  error  propagator  function  can  be  generalized  in  the 
following  way: 
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Error  Propagator  =  Trace 


(*'C7'»  r‘ 


(B-23) 


where  Ct  is  the  correlation  matrix  and  the  A  matrix  is  designed  to  produce  a  full-aperture  lilt  of  zero.  The 
correlation  between  orthogonal  tilts  is  negligible,  but  for  parallel  tilts  between  section*  of  diameter.  D 
separated  by  a  distance,  d,  the  correlation  is  [32] 


Ct(d)  =  -4.873  cos(2^)£^f-r 
+  9.746  cos2(^)£tl^r 


n=0 


n\ 


n  +  -f*  n  +  '6 
n  +  5,  n  +  3,  -n  +  ^- 

n  +  in  +  i 


1  n  +  5,  n  +  3,  -  n  +  j 
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(B-24) 


where  Qd  is  the  angle  between  the  tilt  vector  and  the  separation  vector  This  function  is  independent  of 
the  turbulence  distribution. 

The  error  propagator  has  been  computed  for  a  square  matrix  containing  4,  9,  16,  and  25  aperture 
sections,  with  the  results  shown  in  Figure  B-l.  To  a  good  approximation. 


Error  Propagator  =  0. 1  \(DjDs  )5  3  (B-25) 

over  the  range  studied. 


2  3  4  5 

O/D, 


Figure  B-l.  Error  propagator  for  section-tilts  associated  with  a  multiple  beacon  system.  The  computed  values  can 
be  approximated  by  the  function  0. 1 1  ( D/Ds)5/) . 
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APPENDIX  C— RELATIONSHIPS  APPLIED  IN  THE 
DERIVATION  OF  SYSTEM  PERFORMANCE 


C.l  TURBULENCE  AND  WIND  PROFILES  FOR  PROPAGATION  STUDIES 


The  following  two  models  are  frequently  used  by  the  laser-beam-control  community  to  describe  daytime 
turbulence  conditions: 


SLC-Day  turbulence  model 
C2n(ft)  =0 

=  4.008  X  10-13  A"1054 
=  1.300  X  10"15 
=  6.352  X  10-7  ^-2.966 
=  6.209  X  10-16  A-O-6229 


0  m  <  h<  19  m 
19  m  <  R  <  230  m 
230  m  <  R<  850  m 
850  m  <  k  <  7,000  m 
7,000  m  <  R  <  20,000  m; 


(C-l) 


Hufnagel- Valley  turbulence  model  [43] 

C2(R)  =  5.94  x  iO-53  {will)-  R10  exp(-/i/1000)  +  2.7  X  1(H6  exp(-/i/i500) 

+  A  exp(-/i/100)  .  (C-2) 


The  Hufnagel-Valley  model  has  two  free  parameters,  A  and  w.  The  first  is  normally  set  to  1.7  X 
IO-14  m _2/3,  and  the  second,  which  represents  the  average  wind  speed,  can  be  adjusted  to  achieve  the 
desired  high-altitude  shape.  The  most  common  value  for  w  is  21  m/sec,  which  yields  an  expression 
referred  to  as  the  HV-21  model. 

The  temporal  characteristics  of  turbulence  usually  assumed  to  follow  the  form  proposed  by  Bufton 
[44]  are 


v{R)  =  +  30  expj-[(/i  -  9400)/4800]2 j  ^  (C.3) 

where  the  ground  wind  speed  parameter,  v^,  is  usually  set  to  5  m/sec. 

C.2  HARTMANN  SENSOR  PERFORMANCE  EQUATIONS  FOR  AN  n-DETECTOR 
SUBARRAY 

A  Hartmann  sensor  having  an  arbitrarily  large  linear  dynamic  range  can  be  constructed  by  devoting 
an  n  X  n  matrix  of  detectors  to  each  subaperture  subarray.  This  strategy  is  particularly  useful  for  low- 
repetition-rate  compensation  systems  for  which  the  residual  error  at  each  correction  cycle  is  likely  to  be 
large.  For  a  fixed  dynamic  range,  the  linearity  will  improve  as  the  number  of  detectors  is  increased,  but 
a  good  trade-off  is  achieved  when  the  beam  diameter  is  approximately  equal  to  the  detector  dimension 
Under  these  conditions  a  phase  change  of  one  wave  across  the  subaperture  will  move  the  Hartmann  spot 
the  length  of  one  detector.  Therefore, 

Linear  Dynamic  Range  =  2 rc{n  -  1)  (rad  /  subaperture)  .  (C-4) 
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The  requisite  dynamic  range  for  open-loop  turbulence  correction  can  be  obtained  from  the  phase  structure 
function  by  replacing  the  displacement  parameter  r  with  the  subaperture  diameter  d  .  From  Equation  (2-4), 


DAd,)~6.m(dt/r0) 


5/3 


(C-5) 


which  is  approximately  one  wave  for  a  d/r0  value  of  unity.  To  accommodate  peak-to-peak  fluctuations 
that  might  be  as  much  as  a  factor-of-3  stronger,  a  full  range  of  order-3  waves  would  be  needed;  this 
requirement  could  be  satisfied  with  a  4  x  4-element  detector  array. 

Generalizing  from  Equation  (2-16),  the  gradient  expression  for  an  n  X  n  subarray  illuminated  by 
a  Hartmann  spot  matched  to  the  detector  dimension  is  found  to  be 


8  -  2xJ  1  +  '  1 


(d  \2  £(*'  2  i)ri 


vro; 


i=l 


(rad  /  subaperture)  , 


(C  -  6) 


In 

(=1 


where  ri  represents  the  number  of  photoelectrons  collected  by  the  1th  pixel.  The  photon  noise  near  null 
is  identical  to  that  given  in  Equation  (2-17),  but  the  sensor  noise  expression  is 


i=l 


2n{2i-n  - 1) 


GN 
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where  the  substitution  A/  nfi  —  2^,/*  has  been  made.  Applying  the  identity 


£(2  i  -n- 1)2  =  ±n(n  ~  l)(n  +  1)  shows' that 


i=l 


2  _  4n2  n(n  - 1)(«  + 1) 


j(g  N  ) 

\  e  pe) 
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V2 

V  ro  ) 
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(C-8) 


Combining  this  result  with  Equation  (2-17)  yields  the  total  figure  error  due  to  the  wavefront  sensor. 
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C3  THE  POWER  SPECTRUM  OF  THE  FIGURE  COMPONENT  OF  TURBULENCE 

A  method  for  computing  the  phase  power  spectrum  has  been  described  by  Greenwood  and  Fried 
[21],  who  derive  expressions  as  a  function  of  aperture  position  for  both  tilt-included  and  tilt-removed 
cases.  The  equations  are  greatly  simplified  if  one  uses  the  tilt-included  case,  assumes  that  the  phase  is 
measured  over  a  region  of  zero  extent,  and  allows  the  point  of  measurement  to  lie  at  the  center  of  a 
circular  aperture.  The  combination  of  the  first  and  third  of  these  assumptions  yields  a  result  that  is 
essentially  equivalent  to  the  tilt-removed  case  but  with  a  greatly  simplified  mathematical  structure. 

The  power  spectrum  developed  by  Greenwood  and  Fried  is  represented  by  the  expression 

Ft(f)  =  ~\v/,(T,f)Tt(?)d?  ,  (C  IO) 

where  r  is  the  aperture-position  vector  and  Wp(r,f)  is  the  phase-difference  spectrum.  For  a  circular 
aperture  of  diameter  D, 

vr>(r./)  =  2.079(///0)'83/0-'(r>/'b)53si"-(/V/oD)  •  (C_II) 

where  /0  =  v/nD  represents  a  characteristic  turbulence  frequency  and  v  is  a  wind  velocity  averaged  over 
the  turbulence  path 


v  — 


(C-12) 


The  transfer  function,  7^(r),  comprises  a  set  of  5  aperture-overlap  integrals  that  are  described  in  detail 
in  Greenwood  and  Fried  [21].  For  the  simplifying  assumptions  outlined  above,  only  two  of  these  integrals 
are  non-zero,  which  yields  the  following  representation  of  the  figure  (tilt-removed)  transfer  function: 


cos-'(r/D)-(r/D)[l-(r/D)2 


1/2 


0 


-  n  :  0  <r<  D 
:  D  <  r 


(C  - 13) 


Equation  (C-10)  can  now  be  evaluated  by  combining  Equations  (C- 11)  and  (C-13)  and  employing  the 
approximations  sin2(x)  ®  x 2-x4/3  for  small  values  of  x  =  frlf^p  and  sin2(x)  =  1/2  for  x  large.  The  result  is 
an  integral  over  the  scalar  radius  parameter,  r,  having  the  form 
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which  can  be  solved  with  the  aid  of  the  following  set  of  identities: 


(C-15) 


The  final  solution  has  a  relatively  simply  form  in  the  limit  of  very  low  and  very  high  frequencies, 
and  an  approximation  to  the  complete  function  is  obtained  by  joining  these  two  expressions  at  their  point 
of  intersection: 


r,(/)= 


0. 132  sec(£)*p  Da  ^2/s  V;™  f413 
0.O326sec(C)^r5/3/-8/3 


:  /  <  0. 705  />“  Vq3  5y5/3 
:/>  0.705  D-'rf!Sv3£ 


The  units  are  rad2  of  phase  error  per  Hz,  and  the  integral  of  this  function 


(C-16) 


J  Ff(f)df  -  0.06sec (Qkj  D5/3^  =  014(/)>(,j'  ‘ 
0 


(C- 17) 


is  about  5  percent  larger  than  the  expected  result  for  the  piston  and  tilt-removed  variance  averaged  over 
a  circular  aperture. 


C.4  OPTIMIZATION  OF  THE  SENSOR-NOISE  EXPRESSION 

The  noise  expressions  for  the  phase  and  tracking  sensors  can  usually  be  written  in  the  form 

a2  =axm+  bx~n  ,  (C- 18) 

where  x  an  optimization  parameter,  and  the  variables  a  and  b  are  independent  of  x.  Differentiating  with 
respect  to  x  yields  the  solution 


nM<1/m+n) 

maj 


(C- 19) 

By  combining  Equations  (C-18)  and  (C-19),  the  total  error  for  the  optimal  value  of  x  is  computed: 


~(n\m,m+n 

ftT 

-n/m+n 

\  m) 

Km) 

(C-20) 
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The  relative  contributions  of  the  two  terms  is  as  follows: 


Therefore, 


and 


when  x  has  been  optimized. 


(C-21) 


(C-22) 
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The  application  of  adaptive  optics  technology  to  laser-beam  control  has  been  under  investigation  within  the  military 
community  since  the  early  1970s.  This  work  has  encompassed  a  vide  range  of  theoretical  investigations  as  well  as  an 
impressive  suite  of  laboratory  and  field  experiments.  During  this  period  a  number  of  real-time  systems  incorporating  several 
hundred  degrees  of  freedom  have  been  built  and  tested,  and  the  technology  to  extend  this  number  by  at  least  an  order  of 
magnitude  iB  now  under  development. 

To  exploit  adaptive-optics  concepts  in  the  construction  of  astronomical  telescopes,  the  budgetary  restrictions  under  which 
this  field  of  research  operates  must  be  considered  and  improvements  in  system  performance  must  be  carefully  balanced 
against  equipment  complexity,  reliability,  and  cost.  To  accomplish  this  design  trade-off,  a  set  of  concise  analytic  models 
describing  the  essential  properties  of  a  laser  guide-star  phaBe-conjugation  system  has  been  assembled.  With  the  aid  of  these 
models,  an  optimization  strategy  for  mating  adaptive  optics  to  a  4-m-class  optical  telescope  has  been  developed;  research 
shows  that  such  a  system  might  be  expected  to  improve  the  effective  atmospheric  seeing  conditions  by  nearly  a  factor  of  10 
within  the  isoplanatic  patch  of  the  turbulence  probe. 
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